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ABSTRACT 

We construct the Numerical Galaxy Catalog (lyGC), based on a semi-analytic model of galaxy formation 
combined with high-resolution A^-body simulations in a A-dominated flat cold dark matter (ACDM) cosmo- 
logical model. The model includes several essential ingredients for galaxy formation, such as merging histories 
of dark halos directly taken from iV-body simulations, radiative gas cooling, star formation, heating by super- 
nova explosions (supernova feedback), mergers of galaxies, population synthesis, and extinction by internal 
dust and intervening H I clouds. As the first paper in a series using this model, we focus on basic photometric, 
structural and kinematical properties of galaxies at present and high redshifts. Two sets of model parameters 
are examined, strong and weak supernova feedback models, which are in good agreement with observational 
luminosity functions of local galaxies in a range of observational uncertainty. Both models agree well with 
many observations such as cold gas mass-to-stellar luminosity ratios of spiral galaxies, H I mass functions, 
galaxy sizes, faint galaxy number counts and photometric redshift distributions in optical pass-bands, isopho- 
tal angular sizes, and cosmic star formation rates. In particular, the strong supernova feedback model is in 
much better agreement with near-infrared (/f -band) faint galaxy number counts and redshift distribution than 
the weak feedback model and our previous semi-analytic models based on the extended Press-Schechter for- 
malism. Observed Tully-Fisher relations for bright galaxies and color-magnitude relations for dwarf elliptical 
galaxies in clusters of galaxies are broadly reproduced, but no agreement with observations is obtained over 
whole ranges. Nevertheless, they are improved compared with results produced by other semi-analytic models. 
A direction for overcoming the remaining problems is discussed. We also find that the resolution of iV-body 
simulations, which is down to 3 x 10^ Af© for the minimum mass of dark halos consisting of ten dark matter 
particles in our model, plays a significant role in galaxy formation, and that merging histories of dark halos 
directly taken from A^-body simulations produce results different from models based on the extended Press- 
Schechter model even if the mass function of dark halos at present is set to be the same as that obtained by the 
same A^-body simulations used here. 

Subject headings: cosmology: theory - galaxies: evolution - galaxies: formation - large-scale structure of 



1. INTRODUCTION 

Galaxies are one of the most important hierarchies com- 
posing the Universe. For example, in order to understand the 
evolution of the Universe and its structure, galaxies are often 
used as a measure of space-time in performing wide-field or 
deep surveys. In the field of observational cosmology, there- 
fore, galaxies are indispensable units in the Universe. What 
we should note about galaxies is, however, that galaxies have 
also been evolving along with the evolution of the Universe. 
Thus, correctly modeling the evolution of galaxies is essen- 
tial in observational cosmology. At the same time, galaxies 
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themselves are also interesting objects because many physical 
processes are involved in galaxy formation. However, since 
some processes such as star formation are still poorly under- 
stood, it has been challenging to understand galaxy formation. 
To attack this subject, correctly modeling galaxy formation is 
critical as well, particularly under the current situation that 
numerical hydrodynamical simulation s of galaxy formation 
still h ave large uncertainties (see, e.g.. lOkamoto et aT]l2003l 
l2005h . Thus, the purpose of this paper is to construct a better 
and more reliable model of galaxy formation aided by high 
resolution A-body simulations, which are well established in 
contrast to hydrodynamical simulations of galaxy formation, 
with quality high enough to be compared with recent high 
precision observations. 

It has been widely understood that the cold dark matter 
(CDM) model well describes the Universe. Recent high res- 
olution observations such as those measuring anisotropics 
of the cosmic microwave background by the Wilkinson Mi- 
crowave Anisotropy Probe (Snerael et al. 2003) a nd the type 
la supernova (SN) rate (e.g. iRiess et aljl998tiPerlmutter et alJ 
1999) have revealed that the energy density in the Universe 
is dominated by dark energy or the cosmological constant, 
and that gravity is dominated by the CDM. The large-scale 
structure of the universe has evolved from small density fluc- 
tuations existing in the early universe via gravitational in- 
stability. Because the predicted amplitude of density fluc- 
tuations is monotonically increasing toward smaller scales. 
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small-mass dark halos, which are virialized objects of dark 
matter, firstly collapse and then those of larger masses subse- 
quently collapse swallowing the smaller dark halos. There- 
fore, the formation of cosmological structure predicted by 
the CDM model is often called the hierarchical clustering. 
Most baryonic processes, such as gas cooling and star forma- 
tion, take place within dark halos because of the gravitational 
potential of dark matter. Therefore, galaxy formation mod- 
els must be constructed within the framework of hierarchical 
clustering as a natural consequence of the CDM model. 

Along this line, semi-analytic (SA) model- 
ing of galaxy formation ha s b een d evel oped (e.g. 
Kauffmann. White & Guiderd onil fT993t iCole et alJ 
1994t 120001: Iva n Kampen' Himenez & Peacockl Il999t 
Somerville & Primack 1999; Nagashima, Gouda & Suaiura] 
1999; Nagashima et al. 2001, 2002; Menci et al. 2004 
Hatton et alJ l2003l) . In most SA models, merging histo- 
ries or merger trees of dark halos are realized by using a 
Monte Carlo method based on extended Press-Schechter 
(EPS) models (Peacock & Heavens 1990; Bond et al. 1991'; 
iBower 1991; Lacev & Cole 1993; Kauffmann & White 1993; 
ISomerville & Kolatt 1999; Cole et al. 2000) which provide 
mass functions of progenitor halos. The EPS formalism 
is an extension of the Pre ss-Schechter (PS) model of mass 
functions of dark halos JPress & SchechteJ [l^3)> which 
are derived from the power spectrum of initial density 
fluctuations combined with a spherically symmetric collapse 
model to describe non-lin ear evolutio n of density fluctuations 
JTomital (1969: Gunn & Got3 Il972h . Considering density 
fluctuations within an overdense region which collapses later, 
mass functions of progenitor halos of a larger, later collapsing 
halo are obtained. SA models also include several important 
physical processes such as gas cooling, star formation, 
heating by SN explosions (SN feedback) and galaxy mergers, 
which take place within dark halos. Moreover, to directly 
compare theoretical results with observations, a stellar popu- 
lation synthesis technique is usually incorporated. Because 
we do not follow the dynamics related with these processes 
as hydrodynamical simulations, the computation is much 
faster This means that we can explore a much larger space 
of parameters characterizing these processes, and therefore 
SA modeling is an important approach, complementary to 
hydrodynamical simulations in understanding the physics of 
galaxy formation. 

In addition to those galaxy formation processes, obser- 
vations of high redshift galaxies suffer from several effects 
that are not intr i nsic i n galaxies. Baughetal. (1998) and 
iNayash ima et al.! (1200 U 12002 > considered absorption of emit- 
ted light from high redshift gal axies by intervening H I clouds 
(lYoshii & Peterson>il994i.Madaul 1995ft . The latter work also 
took into account selection effects arising from cosmologi- 
cal dimming of surface brightness dependent on observations 
jYoshii 1993; Totani & Yoshii 2000; Totani et al..2001 ). 

Furthermore, a new ingredient related to SN feed- 
back has recently been introduced into an SA model by 
iNagashima & Yoshii (2004, hereafter NY04) for the first 
time. They considered dynamical response to starburst- 
induced gas remov al from galaxies o n the size and veloc- 
ity dispersion (e.g. lYoshii & Arimotoi ri987). followed by a 
change in the depth of the gravitational potential well. They 
also showed that the model can successfully explain many as- 
pects of photometric, structural and kinematical properties of 
dwarf ellipticals. This is called the Mitaka model, and is the 
basis of the SA model we shall show in this paper. 



The PS model, however, has been recognized to be con- 
siste nt with TV-body simulations only within a factor of 
two (lY ano. Nagashima & Gouda" 199(?; Nagashima & Goudaj 
[ 19971 Monaco 1998; Nagashima 2001; White 200 J 
iReed et alJ 2003). Some SA models, therefore, use 
corrected mass functions based on high resolution iV- 
body simulations (Sheth & Tormen 1999; Jenkins et al. 20011 
'Yahagi, Nagashima & Yos hii,»2004.) instead of the PS model. 
Benson et al. (2000, 2001) used information on halo distribu- 
tion from an A^-body simulation and adapted their SA mod- 
els, GALFORM (Cole et al. 2000), in which an EPS model is 
used. This method has the virtue, not only of allowing spa- 
tial distribution of galaxies to be followed, but also of resolu- 
tion to low-mass halos at high redshift being freely set from 
the A^-body simulations. Unfortunately, we cannot avoid the 
limitations of the PS model with these approaches, because 
the problems of EPS formalism still remain. To bypass the 
limitations of EPS formalism, hybrid SA models have been 
developed in which merger trees of dark halos directly taken 
from iV-body simulations are used, while we have to carefully 
interpret any results from such hybrid models because of the 
limited numerical resolution. Fortunately, currently available 
iV-body simulations are of high enough resolution to resolve 
the effective Jeans mass after the cosmic reionization. This 
kind of hybrid model can also be used to derive spatial distri- 
butions of galaxies. Thus, combined with photometric proper- 
ties, mock catalogs can be made for comparison with individ- 
ual galaxy surveys, without any loss of SA model advantages. 
Therefore, this approach provides strong tools with which to 
study galaxy formation. 

The pioneering work for this approach was carried 
out by Roukemaetal. (1997). Since t hen, their work 
has been substantially improved upon b y iKauffman^^lJ 
( 1999a b), Diaferio et al. (119991 120011). IS omervifl e et alJ 
( 120011) . IHellv et ab (i2003a ^. and 'Hatt on et ah (,2003). We 
have learned many things about the spatial clustering of galax- 
ies, particularly of high redshift galaxies, from those models. 
However, as we shall show in this paper, those models suffer 
from limited numerical resolution of iV-body simulations. As 
another direction to the extension of SA models with A ^-body 
simulations, we should note works by Okamoto & Naga shima! 
( 2001 , 2003) and Sori ngel et al.l ( l2001l) . in which merging his- 
tories of subhalos in galaxy clusters are taken into account. 
This enables us to derive relationships between galaxy proper- 
ties and spatial distribution in halos, such as the morphology- 
density relation (e.g. Dressier 1980). In this paper, we do 
not follow the merging histories of subhalos, but it would be 
worth including those in future. 

In this paper, based on our previous work on the Mitaka S A 
model (NY04) and A^-body simulations using a parallel adap- 
tive mesh refinement (AMR) code (Yahagi 2005), we con- 
struct the Numerical Galaxy Catalog {vGC), which is an SA 
model combined with high resolution A^-body simulations. 
As shown below, the mass resolution is high enough to re- 
solve the Jeans mass, and the size of the simulation box is 
also large enough to e xpress a fair sample of the universe . 
We refer the reader to Yahagi. Nagashima & Yoshii ( 2004*), 
in which mass functions of dark halos are presented, to see 
how powerful this AMR A^-body code is. As a first paper 
in a series of i^GC, here we focus on photometric, structural 
and kinematical properties of galaxies at present and high red- 
shift such as luminosity functions, cold gas fractions, sizes, 
the Tully-Fisher relation, faint galaxy number counts, redshift 
distributions, cosmic star formation histories, and so on. Spa- 



Numerical Galaxy Catalog I. 



3 



TABLE 1 
COSMOLOGICAL PARAMETERS 



0.3 0.7 0.7 0.9 0.048 

Note. — Cosmological parame- 
ters are based on the WMAP results 
ISpergel et alllool . 



tial distribution o f gala xies is discussed in a subsequent pa- 
per (|Yaha2i et aLll2005l). We also plan to extend this model 
in future (Jinoki et al.'vOOS') to include a quasar formation 
model deve loped by Enoki, Nagashima & Gouda (2003) and 
lEnokiet al. (2004). 

This paper is outlined as follows: In §2 we describe our 
SA model. In §3 we summarize model parameters in our SA 
model. In §4 we compare the theoretical predictions of SA 
model galaxies with various observations at present. In §5 
we examine high redshift galaxies predicted by our S A model 
with observations. In §6 we provide summary and conclu- 
sions. 

2. MODEL 

The galaxy formation scenario on which our model is based 
is as follows: In the CDM universe, small dark matter halos 
compared with typical present mass of halos emerge through 
non-linear gravitational evolution of small density fluctua- 
tions. Such dark halos cluster and merge into larger halos 
in a manner that depends on the adopted power spectrum of 
the initial density fluctuations. In each of the merged dark ha- 
los, radiative cooling of virialized hot gas, star formation, and 
gas reheating by supernovae occur. The cooled dense gas and 
stars constitute galaxies. Dark halos are filled with the viri- 
alized hot gas, so if the mass of halos is sufficiently large, it 
should be recognized as an intracluster medium. These galax- 
ies sometimes merge in a common dark halo, and then more 
massive galaxies form. Repeating these processes, galaxies 
form and evolve to the present epoch. 

Merging histories of dark halos are directly taken from iV- 
body simulations ( Yahagi, Nagashima & Yoshii 2004). The 
galaxy formation model running on the merging histories of 
dark halos is almost the same as the Mitaka model given by 
NY04, in which the merging histories were realized by using 
a Monte Carlo method based on the EPS formalism provided 
by Somerville & Kolatt ( 1999). 

Some ingredients of our SA model are revised from those 
in NY04. Modifications include a model of the star formation 
(SF) time-scale, a prescription for mergers between satellite 
galaxies, and a dust extinction model during starbursts. The 
details are described below. 

Throughout this paper, a recent standard ACDM model 
is considered, that is, Hq — 0.3, 57a — 0.7, h — 0.7 and 
(78 — 0.9, where these parameters denote the mean density of 
the universe, the cosmological constant, the Hubble parameter 
and the normalization of the power spectrum of the initial den- 
sit y fluctuatio n field. The shape of the power s pectrum given 
bv ISugivamal iT995ii is used, an extension of iBardeen et al.l 
taking into account the effects of baryons. The cosmic 
baryon density fib — 0.048 is adopted (Spergel et al. 2003|. 
Those parameters are tabulated in Table [2 Model parameters 
described below are summarized in §3. 
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Fig. 1 . — Mass functions of dark halos at ^ = 0. The open circles and solid 
squares with error bars represent inass functions given by the siinulations of 
L = 70h^^ Mpc and 140 h^^ Mpc, respectively. The eiTor bars indicate 1 
cr scatter The minimum masses of those simulations are indicated by Afyo 
and J\/i4o in the figure. T he solid line is a fitting function to results of several 
simulations i Yahaai. Naga shima & Yoshii200'l . 



2.1. N-body Simulations and Construction of Merging 
Histories of Dark Halos 

Merger trees of dark halos are directly taken from cosmo- 
logical iV-body simulations. Our simulation has N — 512'^ 
dark matter particles in a box size of L — 70ft-~^Mpc. Since 
the minimum number of particles identifying a dark halo is 
10, the minimum mass of dark halos il/70 is nearly equal 
to 3.04 x_lO^M0. The simulation code is developed by 
lYahagil (l2?K)l . based on AMR technique to compute the 
gravitational for ce ( Kravtsov et al. 1997; Knebe et al. 20011 
lYahagi & Yoshiili2001.: ■Tevssier..20 02) with periodic bound- 
ary conditions. The method of identifying d ark halos is 
the friends-of-friends (FoF) grouping algorithm dPavis et alJ 
1985), with the linking length h ~ 0.2 as usual. 

Rigorously speaking, however, fully resolved mass might 
be larger, about 5 x IG^^A/q, at which the resultant mass 
function deviates from a power law, as indicated in Figure 
^by open circles with error bars. The solid line in the figure 
is a fitting function for several simulations performed with 
different b ox sizes from 35 to 140 Mpc proposed by 
lYahagi. Nagashima & YoshiH (l2004h . To evaluate the resolu- 
tion effects on galaxy formation, we shall use another set of 
TV-body results of a 140 fe"^ Mp c box with the same number 
of particles, — 512'^ in §Sj2] The mass function derived 
by this simulation is indicated in the figure by solid squares, 
with a minimum mass of Muq — SAfyp. Parameters of those 
two simulations, "70" and "140", are tabulated in Table|2] 

Figure|2]shows the importance of examining the mass reso- 
lution of A^-body simulations on galaxy formation. The solid 
lines indicate masses of dark halos for constant circular ve- 
locities, Kirc ==20, 30, 40, 50, 60, 80, and 100 km s^^ as a 
function of redshift from the bottom to the top. We use the 
spherical collapse model to obtain the relationship between 
the mass and circular velocity, which is expressed as 

M = 3.26x 106A(z)-i/2 ( -Y^\\-^Mq, (1) 

\ km s ^ / 
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TABLE 2 

Parameters of Af-BODY simulations 



Simulation 


L ih^^Mpc) 


N 


Minimum Halo Mass (Mq) 


b 


70 


70 


512^ 


3.04 X l(f(Mro) 


0.2 


140 


140 


512» 


2.43 X 10iO(A/i4o) 


0.2 



Note. — L is the size of the simulation box, A'^ is the number of dark matter particles, 
and b is the linking length parameter for identifying isolated dark halos using the FOF 
algorithm. Particles within b times a mean separation of particles are grouped. The 
minimum halos are defined as those of 10 particles. 
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Fig. 2. — Relationship between mass and circular velocity of dark halos 
as a function of redshift. The solid lines indicate the mass of halos with 
constant circular velocities, Vcirc = 20, ...100 km s~^, from the bottom to 
the top, as indicated in the figure. The dashed lines represent the minimum 
masses A/70 and A/140 of the simulations with 70 and 140/i~^ Mpc box 
sizes, respectively. 



where A(2) is the ratio of the mean density in halos collapsing 
at z to the critical density of the universe. The dashed lines 
are the minimum masses Afyo and M140 as indicated in the 
figure. 

After the cosmic reionization, the Jeans mass should be at 
least T4irc — 20 km s^^, which correspon ds to T 10^ 
K. Based on his hydrodynamical simulations, 'Gnedin' C200Cf) 
suggested that an effective Jeans mass, or the filtering mass, 
becomes even larger affected by the past thermal history, al- 
though the results for low redshifts were an extrapolation 
aided by a linear perturbation theory. Somerville (2002) has 
reported the following result on the effective Jeans mass: the 
mass of halos which retain half of the gas is well approxi- 
mated by mass corresponding to circular velocity Vcirc — 50 
km s^^; halos with VJ;iic < 30 km s^^ are not able to ac- 
crete most of the gas and those with V^iic > 75 km s^^ are 
almost unaffected. As shown in Figure|2j Afyo corresponds to 
Vcirc — 20 km s"^ at z ~ and 40 km s^^ at 2: ~ 5. These 
values should be sufficient to resolve the effective Jeans mass. 
On the other hand, il/140 corresponds to Vcirc — 40 km s^^ 
at z ~ and 80 km s^^ at z ~ 5. These are marginal, unless 
the effective Jeans mass is similar to or larger than that es- 
timated by Gnedin (2000). Thus our 70 /i^^Mpc simulation 
is quite plausible for correct modeling of galaxy formation. 
Apart from detailed discussion on the Jeans mass, we shall 
study the effects of the resolution on luminosity functions in 



§Sj2] Note that other theoretical data in the public domain 
provided by SA models combined with A^-body simulations 
have much larger masses of dark matter particles. For ex- 
ample, the particle mass in the GIF project (Kauffmann et aQ 
'1999a) is about 65 times larger than ours, and that of the 
GallCS project (Hatton et al. 2003) is about 27 times larger 
These masses are even larger than our 140 h^^ Mpc box sim- 
ulation. 

In our simulation, first collapsing halos emerge at z ~ 18. 
Then we stock data sets of the particle distribution for ap- 
proximately 70 time slices. When identifying dark halos, a 
marker particle for each halo is also defined, which is a par- 
ticle whose mechanical energy is minimum. The marker par- 
ticle is assigned to the halo's central galaxy, which is defined 
below. The position of galaxies is specified using the corre- 
sponding marker particles. A halo at time ti is connected with 
another halo at the next time-step fj+i (> ti) as its progenitor 
halo, which has the largest number of dark matter particles 
common to the both. 

Some collapsed halos, especially small halos, disappear 
later. Such "floating" halos should be fake due to fluctuations 
of particle distribution, so we remove the floating halos from 
the sets of merger trees. Some collapsed halos sometimes 
fragment, particularly at low redshift. This is caused by highly 
dense clumps penetrating their host halos, which go back into 
and are absorbed by the host halos again later. To avoid such 
a "loop" structure in merging histories, the fragmented halos 
are forcedly merged into the host halos. Then each merger 
tree has a tree structure with a root halo at z = 0. Almost all 
fragmented halos are small and less than 10% of their host ha- 
los in mass. Because the resolution of our A^-body simulation 
is very high, small halos swallowed by larger halos likely sur- 
vive and penetrate. Physically, gas in such clumps should not 
fragment and should quickly merge into hot gas in their host 
halos, apart from dark matter, because gas is collisional dif- 
ferent from collisionless dark matter Therefore, as far as we 
focus on the evolution of baryons, this manipulation should 
be performed. Th e details of this simulation can be found in 
Yahagi (2005|) and lYahagi. Nagashima & Yoshiil ( ^004) . 

2.2. Tidal Stripping of Subhalos 

Recent high resolution iV-body simulations of galaxy clus- 
ters have revealed that swallowed dark halos survive in their 
host halos as subhalos. Envelopes of these subhalos are 
stripped by tidal force from the host halo. In our simulations, 
such subhalos are not identified but are taken into account 
by the following model. We assume the radius of a tidally 
stripped subhalo rt, which had a radius of and mass of Ms 
before the stripping, by 

^pcri 

rs ■Tapo Ms VKiic,h/ 

where rpdi and Tapo are respectively the pericenter and apoc- 
enter for the orbit of the subhalo, and subscripts "h" and "s" 
indicate the host halo and subhalo. In thi s paper, a ratio of 
''oori/^'aDo ~ 0.2 is assumed, as in NY04 (^ Okamoto & Habel 
1999, 2000). Now we assume a singular isothermal distribu- 
tion for dark matter in subhalos. This means that the decrease 
in mass by the tidal stripping is proportional to rt/rg. Thus, 
the tidal str ippin g significantly affects the time-scale of galaxy 
mergers (§E11'- 

2.3. Hot Halo Gas and Its Cooling 
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Firstly, we define the formation of halos as follows: Dark 
halos grow as time passes via mergers and accretion of other 
halos. Hereafter, the formation epoch of halos is referred to as 
the time when the mass of a dark halo exc eeds twice the ma ss 
of the halo at the final formation epoch (Lacev & Cole 1993|). 
At this time, the circular velocity of the halo is reestimated, 
since the halo is regarded as collapsing at this epoch. The 
mean mass density in dark halos is assumed to be proportional 
to the cosmic mean density at the formation epoch following 
the spherically symmetric collapse model, independent of the 
iV-body results. Each collapsing dark halo contains baryonic 
matter with a mass fraction ili^/flo. The baryonic matter con- 
sists of diffuse hot gas, dense cold gas, and stars. 

When a halo forms with a circular velocity V^iic, the hot gas 
contained in the halo is shock-heated to the virial temperature 
of the halo, 



(3) 



where Wp, /cb and /i are the proton mass, the Boltzmann con- 
stant, and the mean molecular weight. The hot gas density 
is assumed to have an isothermal density profile with a finite 
core radius, according to Shimizu et al. (2002), 



Phot{r) 



Phot,0 



1 + {r/r,y 



(4) 



where rc = 0.22i?vir/c, i?vir is the virial radius of the host 
halo, and c is the concentration modeled by Bullock et al.. 



c{M,z) = 



M 



-0.13 



(5) 



1 + z \1.5x IO13/1-IM0 

While in ICoie et al.l (j2000|), the concentration is altered de- 
pending on the history of each dark halo, we simply assume 
the above form of the concentration. 

A part of the hot gas cools and accretes to the disk of a 
galaxy until subsequent collapse of dark halos containing this 
halo. The amount of the cooled gas is defined by the mass 
of the hot gas encircled by the cooling radius, at which the 
time elapsed from the formation epoch is equal to the cooling 
time-scale derived using the above density distribution. 



t i(r) — '^ Phot(?^) fcer yir 

2 /iTOp n2(r) A(r, 



(6) 



where (r) are the electron density of hot gas at r, Z^^t is the 
metallicity of hot gas, and A is a metallicity-dependent cool- 
ing function provided by Sut herland & Dop ita ( 1993). For 
simplicity, gas which cools in a time-step within a halo is 
added to col d ga s of the central galaxy of the halo, which is 
defined in i! ^2.5\ at the beginning of the time-step. Chemical 
enrichment in hot gas is consistently solved with star forma- 
tion and SN feedback, as shown below. 

In order to avoid the formation of unphysically large galax- 
ies, the cooling process is apphed only to halos with Vcirc < 
Vcut- Although iCole et alTi2000l) apparently succeeded in 
avoiding the formation of monster galaxies without such a 
cooling cut-off by introducing a history-dependent core ra- 
dius, they have found that, in the case of higher baryon 
density, f^h — 0. 04, too many monster galaxies emerge 
jBenson et alJ2003h . Here we take a simple approach to avoid 
the formation of monster galaxies, which is stopping cooling 
by hand. This will be solved if we treat heating processes in 
clusters of galaxies correctly, e.g. heating by AGN/QSOs. 
Throughout this paper, we set Vcut — 210 km s^^. This 



value might seem somewhat small. In addition to the higher 
baryon density, the reason for the difference of this kind of 
cooling cut-off from ones used in other recent SA models, in 
which Vcut ^ 350-500 km s^^ are used, is the treatment of 
reheated gas by SNe. For example, superwinds are adopted in 
iBaugh et al . (2005) and Na gashima et al.. (J005a) which ex- 
pel gas by SNe in addition to the usual SN feedback. This 
frees us from the artificial cooling cut-off. However, the same 
model, but without the superwind, must ad opt a small value 
of circular velocity for the cooling cut-off. Nag ashima et alJ 
( 2005a ) have reported that a model without the superwind re- 
quires V^ut = 100(1 + z)^/^ km s^^, where the redshift de- 
pendence originates by assuming heating by conductive heat 
fluxes from outer envelopes of halos. Thus, in the case of 
r^b — 0.048, under the assumption that all baryons are con- 
tained in halos, such a small Vcut is required. In future, mech- 
anisms of gas-heating and of gas-expulsion from halos should 
be modeled. 

The above assumption that all baryons are retained by dark 
halos might be too simple. Some have argued that gas known 
as the warm/hot intergalactic medium (WHIM) may be out- 
side halos (e.g.lC en & Ostriker illl99^lYoshikawa et alJl20Wt 
iNicastro et al.l l2005). If the WHIM is the gas expelled from 
galactic-scale halos caused by the superwind or similar mech- 
anisms, the baryon fraction in the halos is lowered. In- 
stead, the WHIM might be the gas contained in smaller ha- 
los which are not re solved in hydrodynamical simulations 
(e.g. llnoue et al.l2004 ). Future observation of the WHIM will 
provide a clue to the baryon fraction in galactic-scale halos, 
which constrains the value of Vcut- 

2.4. Star Formation and Supernova Feedback 

Stars in disks are formed from the cold gas. The SF 
rate (SFR) ■0 is given by the cold gas mass Mcoid and a 
SF time-scale as ?/; = Mcoid/T"*- Although determining 
T^: in a purely theoretical way is quite difficult, recent SA 
models suggest that should be dependent on character- 
istic velocities of galaxies and independent of redshift to 
explain many observations ( Kauffmann & Haehnelt 200(1 
Nagashima et al. 2001, 2002; Enoki, Nagashima & Gouda 
2003; Okoshi et al. 2004; Nagashima & Yoshii 2004 
Baugh etalJ l2005t iNagashima et alJ l2005ah . In this paper 
we adopt a simple prescription, similar to that of ICole et alJ 
ri994) . 

,fVa-"' 



(7) 



where and a, are free parameters, chosen to match the 
ratio of the cold gas mass to the _B-band luminosity of spiral 
galaxies and Vd is the disk rotation velocity defined in Wl.6\ 
Vhot is a parameter related to the SN feedback (see below). In 
NY04, the star formation time-scale was 



Vhot 



(8) 



where ahot is also a parameter related to the SN feedback 
defined below. Because this prescription has only one free 
parameter, r^, it is very helpful to search for suitable param- 
eters. Although this worked well in NY04, in fact, we found 
that this cannot reproduce cold gas fractions in spiral galaxies 
in the case of weaker SN feedback. 

Massive stars explode as Type II SNe and heat up the sur- 
rounding cold gas. This SN feedback reheats the cold gas at 
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a rate of Mrohcat = Mcoid/Trohcat, where the time-scale of 
reheating is given by 

'^reheat — n/T r \ ' 



where 



VHot 



(10) 



The free parameters Vhot and ahot are determined by match- 
ing the local luminosity function of galaxies with observa- 
tions. 

Below, we shall use two parameter sets. One of them has 
a strong SN feedback efficiency with ahot = 4. This model 
is hereafter called the "strong feedback" (SFB) model. In this 
case, because of the strong SN feedback, we can take a* = 
— tthot, which makes a star formation time-scale similar to 
NY04. It is not the case for a weak SN feedback. In the other 
model, the SN feedback is weak, ahot — 2. This is called the 
"weak feedback" (WFB) model. 

With the above equations and parameters, we obtain the 
masses of hot gas, cold gas, and disk stars as a function of 
time or redshift. 

Associated with star formation and SN feedback, the 
che mical enrichrne nt is also taken into account by extend- 
ing "Maeder ( 1992). For simplicity, instantaneous recycling 
is assumed for SNe II, and any contribution from SNe la 
is neglected. The recent progress on the chemical en- 
richment with SNe la in the framework of SA models is 
found in Naaashima & Okamoto ( 2004) and Naga shima et al... 
(|2005ab). 

To sum up the above, the basic equations are 

M,=aV', (11) 
Mcoid = -(a + /3)V', (12) 
Mhot=/3V', (13) 

^cold 

(14) 

(Mhot^hot)=/3^cold^, (15) 

where the dot ( ) indicates the time derivative, Af, and 
Mhot are the masses of stars and hot gas, respectively, ip = 
-^com/t* is the SFR, a is a locked-up mass fraction which 
is set to be 0.75, consistent with a stellar evolution model 
we use, Zcoid and Zhot are the metallicities of cold and hot 
gases, respectively, and p is the chemical yield which is set 
to be twice the solar metallicity in order to be consistent with 
the observed metallicity distribution of elliptical galaxies. We 
can solve these equations analytically as 



A A/cold 
AAf,(i): 

AMhotW^ 

Zcold{t) - 
Zhot {t) - 



:Af" 



)ld 



a 



a + i3 
f3 



a 



f3 



1 — exp 
AMcoid(0, 
AMcoid(0, 



-(a 



(16) 



- Zco\d 



t 



P- 



(18) 



(19) 



-'"hot'^hot 



-f3 



{ 



Zcoid(i) A A/cold (i) 



a + (3 

{Z,Mt) - ^coid)Mcoid }]/Mhot(i)(20) 



where variables with A indicate increments or decrements 
which are defined as positive and the superscript stands for 
initial values at the beginning of the time-step, at which time 
t in the solutions is set at zero. For reference, we describe a 
mass-weighted mean stellar metallicity because it should be 
helpful to see the role of SN feedback on the chemical enrich- 
ment. The mass-weighted metalhcity is given by 



/* AUZcoid{t)dt 



Z 



P 



1 



ue 



old 



(21) 



a + j3 1-e-" 
where u= {a + 0)1/ t^. This means that when the SN feed- 
back is strong, the increase in mean stellar metallicity is small. 
The SN feedback clearly plays an important role in chemical 
enrichment. 

2.5. Mergers of galaxies and formation of spheroids 

When two or more progenitor halos have merged, the newly 
formed larger halo should contain at least two or more galax- 
ies which had originally resided in the individual progenitor 
halos. By definition, we identify the central galaxy in the new 
common halo with the central galaxy contained in the most 
massive of the progenitor halos. Other galaxies are regarded 
as satellite galaxies. These satellites merge either by dynami- 
cal friction or random collision. The time-scale of merging by 
dynamic al friction is given by r^j^y^ = /mrgTfric^ where Tfiic is 
given bv lBinnev & Tremaind ((l987), 

1.17 i??Vcirc 

Tfric = , ■ , (22) 

In Ac GMsuh 

where Rh is the radius of the new common halo, Afgub is the 
mass of the tidally truncated subhalo, and Ac is the Coulomb 
logarithm. The parameter /mrg is set to 0.8 in this paper by 
matching the bright end of luminosity functions. In order to 
judge whether or not a satellite merges with a central galaxy, 
we compute the time elapsed for each satellite galaxy since 
the galaxy becomes a satellite. This elapsed time is reset to 
zero if the mass of the host halo doubles. When the elapsed 
time exceeds r,ni g, which is estimated from masses of the host 
halo and a tidally-stripped subhalo harboring the satellite, the 
satellite galaxy is accreted to the central galaxy. 

Satellite galaxies sometimes merge in a time-scale of ran- 
dom collision. The merger rate, k, can be written using phys- 
ical parameters of galaxies and their host halo. 



naVci 



(23) 



where n is the number density of galaxies, and a the cross 
section. Since the cross section is estimated from gravitational 
focusing, we can write cr cx {Grng/V^^j.^)'^, where nig is the 



mass of each galaxy. Using the virial relation, m 



9^1 



/G, 



/ITS 1^ J t= ' a a yi ' 

' ' where Tg is the size of galaxies and Vg a characteristic velocity 



of galaxies, the merger rate is 



k cx 



(24) 



where N is the number of galaxies in the halo. iMakino & Hull 
C1997) derived a merger rate in a situation in which the mass 
of all galaxies was the same, 

N /lMpc\' ( rg 



500 



100km s- 



4 / 300km 



Gyr-\ (25) 
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where Ug and ah are the one-dimensional velocity disper- 
sions. Note that the total merger rate is i?tot = Nkym- 
Unless explicitly mentioned, we use this Makino-Hut merger 
rate fo r satellite-satell ite mergers, a s in|Somerville & Primac^ 
(11999). After that, Mamon (2000) has proposed relaxing the 
assumption of equal-mass mergers in the following way: Con- 
sider a merger of galaxies with masses m and Am and assume 



Because k is considered to scale as 



,,2\2 



the 



rate should be 



fc(m, Am) = k 



l + Ai/3\' /1 + A2/3 



(26) 



To introduce this into our model, the merger rate of a galaxy 
of mass m with galaxies of masses A^m, where i stands for 
other satellite galaxies in the same halo, is 



N 



(fc(m)) = 



fc(m, A^m) 



N 



(27) 



Later we shall see how this effect of non-equal mass merg- 
ers affects our results in Figure [TO] In most of our results, 
however, we do not take into account this effect because it 
is almost negligible, at least at present, and because this re- 
quires quite a long time to compute. Below, we take /cmh 
as a fiducial merger rate, fc, for equal-mass mergers, and use 
quantities of subhalos for the mass, radius and velocity dis- 
persion of satellite galaxies. Here we define a collision time- 
scale, Tcoii = ^/{k)- With a probability A^/tcou, where At 
is the time-step corresponding to the redshift interval Az of 
merger trees of dark halos, a satellite galaxy merges with an- 
other satellite picked out randomly. 

Consider the case when two galaxies of masses mi and 
m2(> mi) merge. If the mass ratio / = mi/m2 is larger 
than a certain critical value of /bulge, we assume that a star- 
burst occurs and that all of the cold gas turns into stars and 
hot gas according to the same SN feedback law, ea.dlOL and 
all of the stars populate the bulge of a new galaxy. Note that 
when appl ying the SN feedback law, Vd is replaced by Vb de- 
fined in §§221 which is the velocity dispersion of the new 
bulge. The time-scale of the starburst is given by the dynami- 
cal time-scale, r» = re/Vb, wh ere is the effective radius 
of the bulge, also defined in ii il2.6l On the other hand, if 
/ < /bulge, no starburst occurs, and a smaller galaxy is sim- 
ply absorbed into the bulge of a larger galaxy. Note that in 
NY04 the smaller galaxy was absorbed into the disk of the 
larger galaxy. We have changed it to resemble the observed 
distribution of bulge-to-disk ratios of galaxies. Throughout 
this paper we use /bulge ~ 0-3, which gives a consistent mor- 
phological fraction in /s u-band galaxy n umber counts [see, 
e.g.. Figure 13 in Nagas hima et al.l ( I2002^ 1. 

2.6. Size of Galaxies and Dynamical Response to 
Starburst-induced Gas Removal 

Here we show a full description to estimate the size of 
galaxies, the rotation velocity of disks, and the velocity dis- 
persion of bulges. The procedure is almost the same as that 
given by NY04 and shown schematically in Figure |3] A new 
parameter providing a fraction of energy dissipation during 
the major merger is introduced. It should be noted that the 
inclusion of the dynamical response to gas removal induced 
by starbursts provides much better agreement with observed 
kinematical and structural parameters, particularly for dwarf 
galaxies, as shown in NY04. 



Dark halos have been considered to acquire angular mo- 
menta during the linear regime of their collapse through 
tidal torques made by the initial density fluctuation field 
JPeebles 1969; White 1984; Catelan & Theuns 1996a bl 
Nagashima & Gouda 1998). Assuming that the hot gas com- 
ponent that cools and contracts within a dark halo has the 
same specific angular momentum as the dark halo has, we 
can estimate a radius at which the gas is supported by rotation 
with the conservation of its specific angular momentum. In 
this way, we determine the sizes of spiral galaxies. 

In this model, angular momenta of dark halos are given by 
the following procedure, independent of the A^-body results: 
This is not only for simplicity, but also because it is challeng- 
ing at this stage to deal with the change in direction of the an- 
gular momentum, which is neglected here as usually assumed. 
The distribution of the dimensionless spin parameter Ah, 
which is defined by Ar = L\E\^/'^ /GM^/'^ where L is the 
angular momentum and E is the binding energy, is well ap- 
proximated by a log-normal distribution ( Mo, M ao & Whit3 
fT998l) . 



p(Ah)c?Ah 



1 



27rCTA 



exp 



(InAH -InA)^ 



d In Ah , 



(28) 

where A is the mean value of the spin parameter and (j\ is its 
logarithmic variance. We adopt A = 0.05 and ax = 0.5 ac- 
cording to Mo. Mao & Wh ite ( 1998). When the angular mo- 
mentum of the hot gas is conserved, the effective radius re 
of a resultant cold gas disk is related to the initial radius Ri 
of the hot gas sphere via r, = (1.68/\/2)AHi?i ("Falll lI979l: 
Fall & Efstathiou 1980; Fafl 1983). The initial radius Ri is 
set to be the smaller of the virial radius of the host halo and 
the cooling radius. A galaxy disk grows due to cooling and 
accretion of hot gas from the more distant envelope of its host 
halo. In our model, the disk size of each central galaxy is es- 
timated at each time-step by the above method, but it is not 
always renewed. Only when both mass and estimated size of 
each disk are larger than those in the previous time-step, the 
size is renewed. At this time, the disk rotation velocity Vd is 
set to be the circular velocity of its host dark halo. 

Size estimation of high-redshift spiral galaxies, however, 
carries uncertainties because of the large di spersion in their 
observed size distribution. For example, Simard et al.l (Il999h 
suggests only mild evolution of disk size against the redshift, 
taking into account the selection effects arising from the de- 
tection threshold of surface brightness, although the above 
simple model predicts disk size proportional to virial radius 
i?vir of host dark halos evolving as i?vir oc 1/(1 + z) for a 
fixed mass. Allowing for the possibility that the conservation 
of the angular momentum is not complete and/or for other 
sources such as the SN feedback to stop the contraction of 
gas disks, we generalize this size estimation by introducing a 
simple redshift dependence, 

1 f)8 

r,^—X^R,{l + zY, (29) 

where p is a free parameter We simply use p = 1 as a refer- 
ence value in this paper in order to be consistent with the ob- 
servational galaxy number counts. The changing of the value 
of p affects the selection effects due to the cosmological dim- 
ming of surface brightness and the dust extinction, because 
the dust column density also changes with galaxy size. 

In the case of the formation of elliptical galaxies driven by 
major mergers of galaxies, the size and velocity dispersion 
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angular momentum and 
central and global velocity dispersions of (sub-)dark halos 

'^Hf ^i) " cent, v " circ 



gas cooling 



disk size and rotation velocity 

re, Vd 



o 
> 

o 

aj 

3 

O 



Other 
galaxies 



major merger, energy dissipation, 
and virialization 



bulge size and velocity dispersion 
before starburst 

ri, Vb,i 



dynamical response to gas removal 



bulge size and velocity dispersion 
after starburst 



dynamical response to gas removal 



central circular velocity of (sub-)dark halos 



V. 



centyf 



disk formation 



— bulge formation 



response to dark halos 



Fig. 3 . — Diagram showing the flow to estimate sizes and velocities of dark halos, disks, and bulges. 



are estimated, explicitly taking into account the dynamical re- 
sponse to starburst-induced gas removal. The outhne is as fol- 
lows: At first, two or more galaxies merge. During the merger, 
a fraction /diss of the total energy of baryonic matter is as- 
sumed to be dissipated and a spheroidal merger remnant is 
formed and reaches the virial equilibrium immediately. Sub- 
sequently, the starburst occurs and a fraction of gas is grad- 
ually expelled from the spheroidal system according to the 
SN feedback law. Because the gravitational potential varies 
during the gas removal, the system expands and its velocity 
dispersion is lowered. The detail is shown below. 

Sizes of early-type galaxies are primarily determined by 
the virial radius of the baryonic component. When a major 
merger of galaxies occurs, assuming some fraction of the en- 
ergy loss to the total energy and no rotation of the merger rem- 
nant, we estimate the velocity dispersion of the merged sys- 
tem. Now we assign the subscript to the merged galaxy, and 
subscripts 1 and 2 to the central and satellite galaxies, respec- 
tively, in the case of central- satellite merger, or to larger and 
smaller galaxies, respectively, in the case of satellite- satelhte 
merger. Using the virial theorem, the total energy for each 
galaxy is 



1 



[MbV^ + {Mi + M,M)V^], 



(30) 



where M(, and are the masses of bulge and disk stars, re- 



spectively, and Vh and Vd are the velocity dispersion of bulge 
and the rotation velocity of disk, respectively. Assuming the 
virial equilibrium, the binding energy E), between the progen- 
itors just before the merger is given by 



Then we obtain 



E1E2 



{M2/Mi)Ei + {Mi/M2)E2 ' 



/diss (-El + E2 + Eb 



(31) 



(32) 



where /diss(> 1) is a parameter which describes how large a 
fraction is dissipated from the system. Just after the merger, 
there remains only the bulge component consisting of cold gas 
and stars in the merger remnant, whose velocity dispersion is 
directly estimated from the above equation. Then the size of 
the system just after the merger is defined by 



GMj 
2^' 



(33) 



where = Af* + Mcoid is the total baryonic mass of the 
merged system. In reahty, the rate of energy dissipation de- 
pends on complicated physical processes such as the escape 
of high velocity stars and the viscosity and friction due to gas. 
Throughoutthispaper, however, we adopt /diss = 1.3 for sim- 
plicity, which is chosen so that the size and velocity dispersion 
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of ellipticals agrees with observations. Apart from /diss, this 
way of estimating Vb and is similar to that adopted in other 
SA models such as that bv lCole et alJ J200Q !). in which they 
formulate the conservation of energy in terms of the size. 

Next, the cold gas turns into stars and hot gas. Since part 
of the gas is expelled from galaxies to halos due to the SN 
feedback, the mass of the system changes. Assuming the den- 
sity distributions of baryonic and dark matters, the dynami- 
cal response on the structural parameters to the mass loss can 
be estimated. This process was taken into account for the 
first time in the Mitaka model (NY04) in the framework of 
SA m odeling. In this paper we adopt the Jaffe model (Jaf^ 
Il983h for baryonic matter and the static singular isothermal 
sphere for dark matter, and assume slow (adiabatic) gas re- 
moval compared with the dynamical time-scale of the system. 
Note that Nagashima & Yoshii (2003) also show the dynam- 
ical response for different density distributions of baryonic 
and dark matter and for rapid gas removal cases. Defining by 
A4,R,Y and U the ratios of mass, size, density and velocity 
dispersion at a final state relative to those at an initial state, the 
response under the above assumption is approximately given 
by 

^-^^/^ (34) 



M + D/2' 



lYR^ + Df{zf)/2 



(35) 



Vb,. V 1 + Df{zi)/2 

where Ai = YR^, D = 1/yizf, y and z are the ratios of 
density and size of baryonic matter to those of dark matter 
and, f{z) is 

ln(l + z) . V 



In 1 



(36) 



The subscripts i and / stand for the initial and final states 
in the mass loss process. The contribution of dark mat- 
ter is estimated from the central circular velocity of halos, 
Vcent, which is defined below. Acco rding to num erical sim- 
ulations of starbursts p erformed by iMori et al.i ( Il997h and 
iMori. Yoshii & Nomotol(IT999) . the time-scale of gas removal 
might be comparable to the dynamical time-scale for dwarf 
ellipticals. If the gas-removal is almost instantaneous, the ef- 
fect of the dynamical response becomes stronger, as shown 
by Nagashima & Yoshii (2003). Therefore, the assumption of 
the slow gas removal is conservative. In reality, some galax- 
ies might become gravitationally unbound due to the rapid 
gas-removal. Stars of such galaxies should contribute to halo 
stars or intracluster stars. 

Finally the back reaction to the gas removal by the SN feed- 
back on dark matter distribution is computed. It is natural to 
consider that dark matter distribution in the central region of 
dark halos is also affected by the dynamical response to gas 
removal, because baryons usually condense with density com- 
parable to that of dark matter. To take into account this pro- 
cess, we define a central circular velocity of dark halo 14ont- 
When a dark halo collapses without any progenitors, Kent is 
set to Vcirc- After that, although the mass of the dark halo 
grows by subsequent accretion and/or mergers, V^cnt remains 
constant or decreases with the dynamical response. When the 
mass is doubled, Vcont is set to Vdrc again. The dynamical 
response to mass loss from a central galaxy of a dark halo by 
SN feedback lowers V^cnt of the dark halo as follows: 

Kent,/ _ M//2 + Md(r,/rd) ^^^^ 



The change of Kent in each time-step is only a few per cent. 
Under these conditions, the approximation of static gravita- 
tional potential of dark matter is valid during starbursts. This 
also applies to subhalos. 

Once a dark halo falls into its host dark halo, it is treated 
as a subhalo. Because we assume that subhalos do not grow 
in mass, the central circular velocity of the subhalos mono- 
tonically decreases. Thus, this affects the dynamical response 
later when mergers between satellite galaxies occur. We as- 
sume that the resultant density distribution remains isothermal 
with Kent at least within the galaxy size. 

The details of th e dyn amical response are shown in 
INag ashima & Yoshiil (120031 |2i)04"). The effect of the dy- 
namical response is the most prominent for dwarf galaxies 
of low circular velocity because of the substantial removal 
of gas due to strong SN feedback (Yoshii & Arimoto 198^ 
Nagashima & Yoshii 2004). If the dynamical response had 
not been taken into account, velocity dispersions of dwarf el- 
lipticals would have been much larger than those of obser- 
vations, determined only by circular velocities of small dark 
halos in which dwarf ellipticals resided. For giant ellipticals, 
on the other hand, the effect of the dynamical response is neg- 
ligible because only a small fraction of gas can be expelled 
due to weak SN feedback. We refer the reader to NY04 for 
detailed discussion of the dynamical response. 

2.7. Photometric Properties and Morphological 
Identification 

The above processes are repeated until the output redshift 
and then the SF history of each galaxy is obtained. For 
the purpose of comparison with observations, we use a stel- 
lar population synthesis approach, from which the luminosi- 
ties and colors of model galaxies are calculated. Given the 
SFR as a function of time or redshift, the absolute lumi- 
nosity and colors of individual galaxies are calculated using 
a population synthesis code by Kodama & Arimoto (1997). 
The stellar metallicity grids in the code cover a range from 
=0.0001 to 0.05. Note that we now define the metalHcity 
as the mass fraction of metals. The initial stellar mass func- 
tion (IMF) that we adopt is the power-law IMF of Salpeter 
form (Salpeter 1955), with lower and upper mass limits of 
O.IMq and QQMq, respectively. In some SA models, it has 
been assumed that there is a substantial fraction of invisible 
stars such as brown dwarfs which have masses smaller than 
0.1 A/0. In this paper, however, we do not assume any contri- 
butions of such invisible stars to the stellar mass according to 
NY04. 

The optical depth of internal dust is consistently estimated 
by our SA model. We make the usual assumption that the 
abundance of dust is proportional to the metallicity of cold 
gas, and the optical depth is proportional to the column den- 
sity of metals. Then the optical depth r is given by 



roc 2 (1 + ^) 7 



(38) 



K„ 



M^l2 + Md{n/rd) 



where r^. is the effective radius of galaxies and 7 is a free 
parameter which should be chosen to predict high redshift 
galaxies being consistent with observations. There are large 
uncertainties in estimating the proporti onahty constan t, but 
we adopt about half the value used by Cole et all ( 12000 ) ac- 
cording to NY04, because it predicts too strong an extinc- 
tion to reproduce galaxy number counts. Wavelength de- 
pendence of optical depth is assume d to be the s ame as 
the Galactic extinction curve given bv ISeatonI (Il979h . Dust 
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distribution is simply assumed to obey the slab dust model 
iPisnev. Davies & Phillippsll989h for disks, according to our 
previous papers. 

In this paper, we introduce dust extinction during starbursts. 
We randomly assign the merger epoch within a time-step un- 
der consideration, according to the time-scale of mergers. The 
SF time-scale of the bursting galaxy is assumed to be the dy- 
namical time-scale, r^/Vb- This time-scale enables us to com- 
pute the amount of metals in cold gas at the output redshift or 
the end of the time-step. Assuming the screen dust model, 
the optical depth is estimated. This procedure is important to 
obtain consistency with observational colors of high redshift 
galaxies. If we do not include this, some galaxies right after 
major mergers become too bl ue to be consistent with obser- 
vations. This is discussed in §fj^3 

The origin of 7 is still unknown. We have introduced 
this parameter just as a phenomenological one. In the WFB 
model, we set 7 = 6, but in the SFB model 7 = 1 to decrease 
the optical depth at high redshift, otherwise too few galax- 
ies were formed to reconcile with observed galaxy number 
counts. One of the reasons this is required is that the dust-to- 
gas ratio may evolve with redshift. At high redshift, the ratio 
might be lower than at low redshift. Another possible reason 
is that the chemical yield we adopt is nearly equal to twice 
the solar value, which might be too large. This large value 
is required to be consistent with observed colors of cluster 
elliptical galaxies, which have been considered to have shal- 
lower IMFs of stars than disk stars, so t hat the ch emical yield 
may be larger than the solar ( Arimoto & Yoshiiirr9 86. 1987|). 
If this is true, the value of the chemical yield must be differ- 
ent depending on the environment, that is, the yield for disk 
stars should be approximately solar, but larger than solar for 
spheroidal components. Because 7 is required to decrease the 
optical depth of dust in disks, by adopting a chemical yield 
similar to the solar value, we will be able to set 7 = 0. As for 
recent progress in investigating the chemical enrichment us- 
ing SA models, we r efer the reader to Naaashima & Okamot^ 
j2004 . lBaugh et alJ ( t2005.) and Nagashima et al. (.200 5a). In 
this paper, we do not enter into any details and maintain 7 
as a phenomenological parameter to reproduce observed high 
redshift galaxies. 

Light emitted from high redshift galaxies is also absorbed 
by intergalactic H I clo uds. We include t his effect in our 
SA mo del, according to lYoshii & Petersoiil (|1994) (see also 
lMadaullT995t iMadau et alJll99 6>). The optical depth for a 
source at redshift Zs at a wavelength A is written as 



T\{Zs) = 



dW 



U _ „-s[{l+z)v]N(W)\ 



dzdw r " / ' 

(39) 

where v — c/A, d^N / dzdW is the number of the Lya lines 
per unit interval redshift per rest-frame Lya equivalent width 
W, N{W) is the hydrogen column density which is related 
to W through the curve of growth for the Lya absorption 
line, and s{v) is the sum of absorption cross sections for the 
lines and continuum of the hydrogen atom dependent on the 
Doppler 6-parameter, which is assumed to be & = 20 km s^^. 
The cross sections for the lines and the continuum are ex- 
pressed as 

.,2 



^nn' } — Jnn' 

nripC 



<,40) 




Fig. 4. — Optical depth as a function of observed wavelengtli for various 
source redsliifts, Zs = 1, 2, ...10 from tlie left to the right. Response func- 
tions of BVRi' z' bands for the Suprime-Cam on the Subaru telescope are 
overlaid by the dotted lines from the left to the right. Note that the vertical 
axis for these is in linear scale. 



and 



QAt: mf,e gn{v) 

{1 - e[vn + {v/c)vn - H}e-"'/''di^ 



,(41) 



where 0{x) is the Heaviside step function, r„' is the effective 
damping constant, g„ {v) is the Gaunt factor, and the other no- 
tations have their usual meanings. We consider Lyman con- 
tinuum absorption {n ~ 1) and the first five lines of Lyman 
series, Lya, ... Lye {n — \ and n' — 2, ...6). The line distri- 
bution can be modeled as 



dzdW 



10.7(1 + z) 



2.37" 



-w/w 
W* 



(42) 



where W* = 0^36A , according to the Model 1 in 
lYoshii & Peterson ( 1994). Combined with response functions 
of the filters we use, we obtain the optical depths for individ- 
ual pass-bands. We simply estimate an optical depth for an 
X-band, tx, as 



exp(-Tx) = 



/ dXRx.xe 
fdXR 



(43) 



where Rx.\ is a response function to the X-band filter The 
difference in magnitudes is given by Am = 2.5tx log e. Fig- 
ure |3 shows the optical depths at Zs = 1, ...10, indicated by 
the solid lines, and response functions of BVRi' z' of the Sub- 
aru telescope by the dotted lines, which are used in fjS] 

We classify galaxies into different morphological types ac- 
cording to the B-band bul ge-to-disk lumi nosity ratio B/ D. In 
this paper, following Simien & de Vaucouleurs I 1986). galax- 
ies with B/D > 1.52, 0.68 < B/D < 1.52, and B/D < 
0.68 are classified as elliptical, lenticular, and spiral galaxies, 
respectivel y. Kauffmann, White & Guiderdoni (1993) and 
Baugh, Cole & Frenk ( 1996) showed that this method of type 
classification well reproduces the observed type mix. 

Luminosity profiles of galaxies are assumed to obey the de 
Vaucouleurs' 1/4 law for E/SO galaxies and the exponential 
law for spirals. These profiles are used to compute the isopho- 
tal magnitude, which is obtained by integrating light from 
the galactic center to a radius at which the surface brightness 
reaches a threshold value. If the threshold is smaller in units 
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TABLE 3 
ASTROPHYSICAL PARAMETERS 



Pai'ameter 



SFB WFB 



Annotation 



Observation 



"hot 

Vhot (kms-i) 
Vcut (km s~^) 
a* 

TO(Gyr) 

PiZQ) 
/bulge 
/mrg 
/diss 

A 

p 

7 



4 
140 

210 
-4 
2.4 
2 

0.3 
0.8 
1.3 
0.03 
0.5 
1 
1 



2 
140 

210 
-2.3 

2 
2 

0.3 
0.8 
1.3 

0.03 
0.5 
1 




supernova feedback-related (§§2.3) 

cooling cut-olf (§§2.3) 

star formation-related (§ i|2.4l 

heavy-element yield (5 ^12.41 

major/minor merger criterion (i; il2.51 
coefficient of dynamic al fric tion timescale 
energy-loss fraction (§ i|2.61 

spin parameter distribution (ii i|2.6l 

redshift dependence of disk size (§ t|2.6t 
redshift dependence of dust optical depth ( 



(§E2 



luminosity functions (FigureQ 

luminosity functions (FigurelTl 

cold gas mass (Figure fTTl 

metallicity distribution (not shown) 
morphological counts (not shown) 
luminosity functions ( Figu relTl 

spheroidal size (Figure fl4l and velocity dispersion (Figure lT^ 
disk size (Figure fT3l 

faint galaxy number counts (Figures rT^M9|and|23l 
faint galaxy number counts ("Fi gures 11811 1 9l and l23 1 



Note. — CKhot ^^nd Vhot 'ire the SN feedback-related parameters (eg. 1101 . and are the SFR-related parameters (eq. 0. /bulge is a critical mass-ratio of merging galaxies 
to distinguish majo r an d minor mergers (^ ^2.51 . /mrg is a parameter to tine-tune the time-scale of mergers due to dynamical friction tS ^2.5l . p is th e ch emical yield as defined by 
^aeder ■ 1 992 ) (eq. 1141 . A and ax characterize the distribution function of the dimensionles s sp in parameter, Ah , with a shape of log-normal type (eq. 1281 . Their values are giv en by 
Imo. Mao & White ' 199^. p is a parameter giving a redshift dependenc e of the disk size (eg. 1291 . /diss represents the energy-loss to the total energy during the major merger teg. 1321 . 
7 is a parameter giving a redshift dependence of dust optical depth teg. 1381 . The last column denotes observational quantities which are used to set the values of the parameters of the 
model. 



of surface brightness, e.g., mag arcsec^^, the contribution to 
the isophotal magnitude comes only from the central region of 
the galaxy. Where the threshold is infinity, the isophotal mag- 
nitude is identical to the total magnitude. Thus, the threshold 
depends on the depth of observations. It is very important 
for high redshift galaxies to estimate the isophotal magnitude 
because of the selection effects for detection. The details of 
how to compute the isophotal magnitude are shown in Yoshii 
(0,993), Totani & Yoshii (2000) and Totani et al. (2001J), and 
in iNagashima et al.i ti200L i2002i) in the framework of SA 
models. 

3. MODEL PARAMETERS 

Here we briefly summarize model parameters. As already 
mentioned, we adopt a standard ACDM model. The cosmo- 
logical parameters are f^o — 0.3, f^A — 0.7, h = 0.7 and 
<Ts = 0.9, and the baryon density parameter is 51b = 0.048, as 
ta bulated i n Table[T| suggested by the first-year WMAP results 
( ISpergel et al. 2003). The shape of the power spect rum o f 
density fluctuations adopted here is that of,Sugivama ( 1995'), 
in which the effects of baryons are taken into account to im- 
prove the power spectrum given bv lBardeen et ah (tl98 6). 

Most of the astrophysical parameters are constrained from 
local observatio ns, according to the procedure discussed in 
INagashima et aL (.200 1. 2002) and NY04. The values adopted 
are slightly different from those in our previous papers. This 
is mainly caused by different merger trees which are directly 
constructed in this paper using A^-body simulations. We shall 
show tha t cor rectly constructing merger trees is very impor- 
tant in 8 ^4.11 In the Mitaka model (NY04), we used a mass 
function w hich fits with results of the A^-body simulations we 
used here ( Yahagi. Nagas hima & Yoshiill2004 >). which means 
that mass functions of dark halos are the same at output red- 
shift, e.g. at z = 0. However, we have found that luminosity 
functions in this model are different from those in NY04, in 
which an extended Press-Schechter model was used to con- 
struct merger trees (see Figure |9}. This shows that merger 
trees of dark halos directly affect properties of galaxies. 

Parameters are determined in the following way: Firstly, 
the SN feedback-related parameters, ahot and V^iot are con- 
strained by matching luminosity functions (see Figure 0. 



Roughly speaking, ahot determines the faint-end slope and 
Vhot the magnitude of L* galaxies. Since observational un- 
certainties in estimating luminosity functions still remain, 
we have a degree of freedom in varying those parameters. 
In the Mitaka model, we have used ahot = 4, a param- 
eter also used in this paper and referred to as the SFB 
model. In many other S A models, howe ver, ahnt, = 2 
is often used (e.g. iKauffmann, White & G uiderdonij Il993t 
iomerville & Primack 1999; Cole et al. 2000). Therefore, we 
provide another parameter set, the WFB model, with ahot — 
2. T4iot is determined for luminosity functions to be consis- 
tent with observations under the constraint of ahot- Second, 
SFR-related parameters, and a,, are determined by match- 
ing cold gas fractions in spiral galaxies, Mhi/^b (see Figure 
II 1> . As mentioned, we have found that a^, = —ahot for the 
SFB model, but a^, = —2.3 for the WFB model, which is 
slightly different from —ahot = —2. Finally, other parame- 
ters such as the mean spin parameter are fine-tuned, but they 
are quite similar to physically motivated values. The values 
of these parameters are tabulated in Tabled 

Figure |5] shows the projected density distribution of dark 
matter taken from the L = 70h^^ Mpc simulation at six dif- 
ferent redshifts from z = to 5. It is evident that dark matter 
is clustering more as time passes. Figure |6l shows the same 
distribution as Figure |5] but for galaxies of Mb < —19.5 
(Mb — 5\ogh < —18.7) in the rest frame. The position of 
galaxies is specified by the marker particles, as mentioned in 
it 42.1\ Sizes and colors of symbols represent the rest frame 
magnitudes and colors of galaxies, respectively. In this plot, 
we do not take into account the selection effects and ab- 
sorption of light by intergalactic H I clouds. At higher red- 
shifts, galaxies become bluer By contrast, clustering prop- 
erties of galaxies seem to be almost independent of redshifts. 
This has already been de scribed as "biased" galaxy forma- 
tion, and IKauffmann et al.i (.1999b.) have shown that the cor- 
relation length at which the two-point correlation function of 
galaxies becomes unity does not significantly evolve as dark 
matter does and that the bias of galaxy distribution to dark 
matter grows toward higher redshift. In a subsequent paper, 
we shall investigate clustering properties of galaxies in detail 
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Fig. 5. — Projected distribution of dark matter in tlie L = 70/i~^ Mpc simulation at six redsliifts: z = (top left), 1 (top middle), 2 (top right), 3 (bottom 
left), 4 (bottom middle) and 5 (bottom right). Brighter regions indicate denser regions. 




p[i<"'Uh] c{i|-'t^] ilBi-iii^l 

Fig. 6. — Same as Figurelslbut for galaxies. Plotted galaxies are Mb < — 19.5 in their rest frames. Larger points represent brighter galaxies. Colors of points 
indicate the rest frame B ~ V color of galaxies. 
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Fig. 7. — Local luminosity functions in tlie (a) B and (b) K bands. The 
solid and dashed lines represent the SFB and WFB models, respectively. 
Symbols with error bars in (a) indicate the obse rvational data from APM 
ILovedav et all 19921 filled squares), ESP i Zucca etalIT997t open triangles), 
Durham/UKST (Ratcliffe et al. 1998, crosses), and 2dF iFollces et al. 1999, 

fn circles). Symbols in (b) indicate the data from 2MASS i Kochanek et aia 
11, open circles), 2dF combined with 2MASS iCole et al. 2001, crosses), 
1 et al. l2003^ open triangles), and .Huang et al. 1.2003,- tilled squares). 
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Fig. 8. — Local luminosity functions. Same as FigurelTl but for theoretical 
models plotted. The solid lines represent the same SFB model as shown in 
FigurelTl indicated by SFB/70 in (a). The dashed lines represent a model with 
the same parameters as the SFB model, but for the simulation box size of 140 
h^^ Mpc, which has an eight times larger minimum mass of dark halos. The 
dot-dashed lines represent the Mitaka model given by NY04, but for the same 
parameters as the SFB model. The minimum mass of dark halos in this model 



is given by Vci 



40 km s~ 



JYahagi et alJ2005l) . 

4. GALAXIES AT PRESENT 
4.1. luminosity functions 

Figure shows B and X-band luminosity functions of 
galaxies. The solid and dashed lines represent the SFB and 
WFB models, respectively. Since the SF time-scale hardly 
affects local luminosity functions, they are determined al- 
most solely by the SN feedback-related parameters, ahot 
and Vhot- Symbols with error bars represent observational 
results from the B-band redshift surveys, such as Auto- 
matic Plate Machine (APM; Lovedav et al. 1992), ESO Slice 
Project (ESP; Zucca etal. 1997), Durham/United Kingdom 
Schmidt Telescope (UKST; RatcUffe et al. 1998) and Two- 
Degree Field (2dF; Fol kes et all 11999). and from the K- 
band redshift surveys given by Two Micron All Sky Survey 
(2MASS; Kochanek et al. 2001 ), 2dF co mbined with 2MASS 
dCole e t al. 2001), Bell et aL (»2003i) and lHuang et"al] j2003h . 
The SFB model agrees well with the luminosity functions of 
2dF and Bell et al. (2003), and the WFB model agrees with 
those of ESP and Huang etal.. C2003,) . 

It is interesting and important to assess the effects of nu- 



merical resolution on galaxies. This corresponds to the lower 
limit of the mass of dark halos. The A'^-body simulation we 
use is of a box size of 70 h^^ Mpc. Since all halos must have 
at least ten dark matter particles, the minimum mass A/70 cor- 
responds to 3.04 X IO^Mq. Here we use another result of an 
iV-body simulation with a 140 h^^ Mpc box and the same 
number of dark matter particles. The minimum mass of halos 
is Afi4o = 2.43 X IQIOMq = SAfyo. The resultant luminos- 
ity functions for the 140 h^^ Mpc box simulation are shown 
in Figure|8]for the SFB model, indicated by the dashed lines 
(SFB/140), with the same astrophysical parameters as those of 
the 70 /i^^ Mpc box simulation {SFB/70; solid lines). It is ev- 
ident that this predicts fewer galaxies and a very shallow slope 
at the faint end. One of the reasons is clearly that the mini- 
mum mass of dark halos becomes larger. This is an effect sim- 
ilar to that resulting from the adoption of stronger SN feed- 
back for low-mass halos so that hot gas cannot cool within 
dark halos below Muvi. This means that if one would like 
to connect an obtained parameter set with physical processes 
such as the SN feedback, the numerical resolution would be 
crucial, while we should be able to construct a predictable 
model regarding those parameters as only phenomenological. 
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Fig. 9. — Local luminosity functions. Same as FigureQ but for theoretical 
models. The solid lines represent the SFB model with the same minimum 
mass A'/i4o as the SFB/140 model. Hot gas within halos between Afro and 
Afi4o do not cool at all. The dashed lines represent the same SFB/140 model 
as shown in Figure|8] The dot-dashed lines represent the Mitaka model with 
the minimum mass Muo . 
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Fig. 10. — Local luminosity functions. Same as FigurelTl but for two added 
variant models. The solid and dot-dashed lines represent the SFB and WFB 
models, respectively, the same as those shown in FigureQ The dashed and 
dotted lines also represent the SFB and WFB models, but including the e ffec t 
of non-equal mass mergers between satellite galaxies described by eg. 1271 . 
Note that luminosity functions of models with and without the effect of non- 
equal mass mergers are almost identical. 



Also plotted by the dot-dashed lines is the Mitaka model, but 
with the same parameters as those of the SFB model, which 
are slightly different from the original ones. In this model, the 
minimum mass of halos corresponds to V^iic = 40 km s^^. 
Part of the reason the faint-end slope is steeper than that of 
the SFB model should be that the minimum mass becomes 
smaller at high redshift when stars actively form. This is, 
however, only part of the reason. This is shown below. From 
this comparison, we would like to stress that we should denote 
attention to the minimum mass of dark halos, or the numer- 
ical resolution, which corresponds physically to the effective 
Jeans mass. 

Further investigation of this problem provides interesting 
features. Figure |9] shows luminosity functions for the SFB 
and Mitaka models with the minimum mass Mhq of halos 
indicated by the solid and dot-dashed lines, respectively, as 
well as the SFB/140 model by the dashed lines. At first, it 
is shown that the SFB/70 model with Mhq does not recover 
the SFB/140 model, except for the cut-off at the faint end, 
Mb — blogh ~ —14. The change in the slope caused by 
the change in minimum mass is quite small. Therefore, the 
minimum mass seems to mainly affect only very faint galax- 
ies. Second, the faint-end slope of the Mitaka model with 
Mi4o also hardly changes except for the cut-off. The bump 



at Mb ~ 5 \og{h) ~ —14 is a characteristic luminosity of the 
faintest galaxies determined by the higher minimum mass. In 
other words, galaxies formed in dark halos of the minimum 
mass without any progenitors should have such luminosities 
under the given strength of the SN feedback. Thus it is diffi- 
cult for those to have fainter luminosities. This figure there- 
fore shows again that the numerical resolution significantly 
affects galaxy formation, and that there is a limitation in the 
extended Press-Schechter formalism which gives us merging 
histories of dark halos. 

We would like to stress here that this is not a convergence 
test. If the merger trees included less massive dark halos, the 
results would change further In this case, the faint end of 
the luminosity function would become steeper because such 
small halos produce very faint galaxies. In our current model, 
although no realistic model of the effective Jeans mass is in- 
cluded, the mass resolution corresponding to Af7o should be 
required for realistic modeling, because adopting Afi4o or 
larger minimum mass corresponds to overestimation of the 
effective Jeans mass. It should be very important to model 
the effective Jeans mass by utilizing higher-resolution simu- 
lations in future. 

Finally, we discuss mergers of satellite galaxies. As shown 
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Fig. 1 1 . — Hydrogen mass relative to B-band luminosity of spiral galax- 
ies. The solid and dashed lines represent the medians for the SFB and WFB 
models, respectively, of the distributions of Mn/Lg. The eiTor bars indi- 
cate 25% and 75% percentile points. Symbols indicate the observational data 
from Huchtmeier & Richter ( 1988, open triangles). Sage 1 1993, crosses), and 
l^m'netl 1 2002, solid squares). Note that the data from Huchtmeier & Richtea 
Tl988l) include only H I. 



in eq.j27ll. our code enables us to take into account the effect 
of non-equal mass mergers between satellite galaxies based 
on Mamon (2000), while this is normally switched off. Figure 
nojshows this effect on luminosity functions. The dashed and 
dotted lines indicate the SFB and WFB models with the effect 
of non-equal mass mergers, while the solid and dot-dashed 
lines indicate the SFB and WFB models without this effect, 
which are the same as those shown in Figure It is easily 
seen that the luminosity functions with and without the effect 
of non-equal mergers are almost identical except for the bright 
end. Therefore, we conclude that, fortunately, we do not need 
to consider the effect of non-equal mass mergers. Thus we 
co nsider only equal mass mergers using the rate of kuH given 
bv lMakino & Hull ^1997^ below. 

4.2. cold gas 

The amount of cold gas, which corresponds to the interstel- 
lar medium, is mainly determined by SN feedback-related and 
SFR-related parameters. The former parameters determine 
the gas fraction expelled from galaxies, and the latter, the gas 
fraction that is converted into stars. Because the SN feedback- 
related parameters are already constrained by matching lumi- 
nosity functions as shown in the previous subsection, we can 
use observed gas fractions in spiral galaxies as constraints on 
the SFR-related parameters, a* and r°. 

Figure shows the ratio of cold gas mass relative to B- 
band luminosity of spiral galaxies as a function of their lumi- 
nosity. We assume here that 75% of the cold gas in the models 
is comprised of hydrogen, i.e., Afn = 0.75Afcoid- The solid 
and dashed lines associated with error bars represent the me- 
dians and 25% and 75% percentiles of the distributions for 
the SFB and WFB models, respectively. We have added 0. 1 
mag offsets to the data of the WFB model for error bars to be 
distinguishable. Symbols indicate the observational data of 
individual galaxi es tak e n from Huchtmeier & Richter tl98^ 
open triangles), ISaed J1993L crosses), and iGarnettI J2002t 
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Fig. 12. — H I mass functions. The solid and dashed lines represent the 
SFB and WFB models, respecti vely. The crosse s and open circles indicate 
the ob servational data from Rosenbe rg & Schneide r i2002) and Zwa an et aEI 
120031) . respectively. Again we set Mm = O.TSMcoid for the theoretical 
models. 



solid squares). Note that the data from'Huchtmeier & RichteJ 
11988) include only H I, and that the fraction of hydrogen 
molecule, H2, usually increases with the luminosity. Al- 
though the medians for both models are in good agreement 
with the observations, the error bars for the theoretical mod- 
els seem to be larger than the distribution of the observational 
data. This might suggest the existence of self-regularization 
mechanisms of star formation to develop small scatters of cold 
gas fraction. Nevertheless, we can say that the models agree 
well with the observational data. 

Adopted values of the parameters for the SF model are tabu- 
lated in Table|3l The SFB model has a* = —4, which means 
T^, cx (3. This is similar to the prescription of the SF time- 
scale in NY04, r* oc (1 + /?). On the other hand, in the WFB 
model, we have to set a* = —2.3 7^ — ahot- Unfortunately, 
the simple prescription of NY04 does not work for weak SN 
feedback models. 

Recent radio surveys have provided good data of H I mass 
functions (HlMFs). Figure shows the theoretical predic- 
tion of HlMFs for the SFB and WFB models, respectively. 
Both models provide almost identical HlMFs, because we 
have chosen the SFR-related parameters to be consistent with 
the observed cold g as fractions. Symbols indicate the ob- 
ser yational data from j Rosenberg & SchneideJ ilOO'A crosses) 
and lZwaan et a open circles). Both models reproduce 

the observed HlMFs quite well. Note that we again assume 
A/hi = 0.75Afeoid. 

4.3. galaxy sizes 

Figure shows the effective disk radii of local spi- 
ral galaxies as a function of their luminosities. The solid 
and dashed lines represent the SFB and WFB models, re- 
spectively. Symbols are the observational data taken from 
'Impev et ai] Jl99 6^. To avoid confusion, we have added 0.1 
mag offsets to the WFB model. Here we would like to briefly 
repeat the procedure to determine the disk size. Dark halos 
form with angular momenta randomly assigned with the log- 
normal distribution of the dimensionless spin parameter. Ah 
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Fig. 13. — Disk size of spiral galaxies. The solid and dashed lines repre- 
sent the averages of disk radii of spiral galaxies in the SFB and WFB mod- 
els, respectively. Error bars on these lines are 1 u scatter in predicted sizes. 
To avoid any confusion, 0.1 mag offsets are imposed on the r esults of the 
WFB model. Crosses indicate the observational data taken from llmpev et alj 
11991 . 



(eq. I28I I. Hot gas contained in a dark halo is assumed to have 
the same specific angular momentum as that of the halo. The 
hot gas cools and contracts until the cooled gas is rotationally 
supported. Thereby the size of disks is determined. Good 
agreement of the theoretical models with the observational 
data shown in this figure means that the standard assumption 
we adopt for determining disk sizes works well when the dis- 
tribution of the dimensionless spin parameter is in a reason- 
able range, estimated from the structure formation theory in 
the CDM universe. 

Correctly estimating disk sizes is very important to compare 
theoretical models with observations because it affects many 
observational properties of galaxies s uch as the surface bright- 
ness and the dust column density. Nagash ima et alJ J2001t 
^02) have discussed such effects caused by the uncertainty 
in estimating disk sizes. 

Figure O shows the effective radii of elliptical galax- 
ies against their _B-band magnitudes. The dots represent 
individual elliptical galaxies in the (a) SFB and (b) WFB 
models. The effective radius is estimated from = 
0.744rh ( iNagashima & Yo shii 2003), where rf, is the three- 
dimensional half-mass radius reas.( l33> and ( I34l l1. To make 
figures clearer, only a quarter of all ellipticals are plotted. 
The solid triangles indicate the observational data taken from 
[Bender. Burstein & Faber ( 1992). The theoretical models 
agree well with the observational data of at least brighter than 
Mb — 51og(ft.) ~ —18. At magnitudes fainter than that, 
the fraction of compact ellipticals in the models seems to be 
smaller than that expected from the observations. This would 
be partly caused by the limited resolution. Another possibility 
is that the SF time-scale at high redshift might be slightly too 
long. As clearly shown by NY04, the size of dwarf ellipticals 
strongly depends on the SF time-scale at high redshift. For ex- 
ample, let us consider the SF time-scale simply proportional 
to the dynamical time-scale. In this case, the SF time-scale is 
shortened for higher redshift, and at z = 2, it becomes nearly 
a quarter of that at z = 0. Noting that the SF time-scale also 
determines the rate of gas consumption, therefore, galaxies 
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Fig. 14. — Sizes of elliptical galaxies. The dots represent effective radii of 
individual elliptical galaxies in (a) the SFB model and (b) the WFB model. 
Only a quarter of all elliptical galaxies is plotted. The filled triangles in - 
dicate the observational data taken from fBender. Burstein & FabeJ 11993) . 
The dashed lines denote constant surface brightnesses, 23.5 and 25.5 mag 
arcsec^^, as indicated in panels. 



are statistically gas-poorer than those in the case of the con- 
stant SF time-scale over all redshifts at high redshift, at which 
galaxies frequently merge. Evidently, because the change in 
size and velocity dispersion due to the dynamical response 
becomes smaller, galaxies with small sizes and large veloc- 
ity dispersions increase. As shown by NY04, some compact 
ellipticals represented in FigureOare closely related to ellip- 
ticals formed via gas-poor mergers. Although we are able to 
make such compact ellipticals by slightly shortening the SF 
time-scale at high redshift with an extra parameter relating to 
the redshift dependence, we would like to leave it constant to 
avoid further complexity of the model. A detailed discussion 
about the relationship between the SF time-scale and the size 
and velocity dispersion of ellipticals can be found in NY04. 

The dashed lines represent constant surface brightnesses of 
23.5 and 25.5 mag arcsec^'^ as indicated in the figure. Be- 
cause the detection of galaxies is limited by the surface bright- 
ness, galaxies with surface brightness lower than the detection 
limit, which depends on observations, cannot be detected. Al- 
though our models apparently predict too many dwarf ellipti- 
cals with radii larger than about \ h^^ kpc, they will not be 
observed unless ultra deep observations are performed. Note 
that th e samples of galaxies from Bender. Burstein & Fabe^ 
lfT992l) may be expected to have sufficiently high surface 
brightness to investigate their structural and kinematical prop- 
erties, because the purpose of the paper is to clarify those 
properties. Thus it should not be a problem that many dwarf 
ellipticals with much larger radii emerge. In fact, the selection 
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Fig. 15. — /-band Tully-Fisher relation of spiral galaxies. The solid and 
dashed lines associated with error bars indicate the medians and 25% and 
75% points of the distributions of 7-band magnitudes for the SFB and WFB 
models, respectively. To avoid confusion, 0.1 mag offsets are imposed on 
the results of the WFB model. The observed TF relations are shown as 
a mean relation give n by Mathewson. Ford & Buchhorn 1 1992, thin dashed 
/mel. iPierce & Tull\l>199X thin dot-dashed line), an AGiovanelli et aLtl99^ 
thin dotted line). 



effects brought about by the surface brightness play a signif- 
icant role for high redshift galaxies. This will be shown in 

m 

4.4. velocity-magnitude relations 

Here we investigate velocity-magnitude relations, that 
is, the Tully-Fisher (TF) relation for spiral galaxies 
(jjullv & Fisher 1977) and the Faber- Jac kson (FB) relation for 
elliptical galaxies (Faber & Jackson 1976). 

Figure ^\ shows the /-band TF relation. The solid and 
dashed lines associated with error bars represent the medians 
and 25% and 75% points of the distributions of /-band mag- 
nitudes for the SFB and WFB models, respectively. While 
gas-rich sp irals tend to be more luminous as indicated by 
ICole et all (i2000 ) and NY04, we show only the results for 
all spiral galaxies. We find that, although not shown, the 
TF relation produced only by gas-rich spirals becomes too 



bright, particularly at Vro 



130 km s and too steep at 



Kot ^ 100 km s^ to reconcile with observations in both 
models. The straight lines indicate the results best fit to 
the ob served TF relations of Mathewson. Ford & Buchhorn 
([1992', dashed), 'Pierce & TuUv. Q992, dot-dashed), and 
iGiovanelli et al. ( 1997, dotted). We assume that the linewidth 
W is simply twice the disk rotation velocity, Viot, as usual. 
For model galaxies, we set Vfot = Vd. 

At faster than V^ot ^ 160 km s"^ [log(yrot/km s"^) ~ 
2.2], both the theoretical models agree well with the obser- 
vational TF relations. Besides, the WFB model is marginally 
consistent with some of the observations down to Vi-ot — 100 
km s^^, while the TF relation of the SFB model deviates and 
has a steeper slope compared with the observations. As shown 
in NY04, the slope is mainly determined by cthot because it 
provides the dependence of the strength of the SN feedback on 
the velocity and therefore affects luminosities, particularly of 
dwarf galaxies. The deviation at small velocities might disap- 
pear if the dynamical response on the rotation velocity to gas 
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Fig. 16. — B-band Faber- Jackson relation of elliptical galaxies. The dots 
represent one-dimensional velocity dispersions of individual elliptical galax- 
ies in (a) the SFB model and (b) the WFB model. Only a quarter of all 
elliptical galaxi es is plotted. The filled triangles indicate the observational 
data taken from lBender. Burstein & Fabeil 119921) . The velocity dispersion of 
the models is simply estimated as oio = 14/\/3- 



removal due to the SN feedback is taken into account, which 
is currently applied only to spheroidal formation by starbursts. 
Because dwarf galaxies significantly suffer from the SN feed- 
back, disks may expand and rotate slowly. 

Figure^]shows the /3-band FJ relation of elliptical galax- 
ies. The dots represent individual ellipticals in the (a) SFB 
and (b) WFB models. The model prediction s broadly agree 
with the observational data from Bender. Burstein & FabeJ 
( 1992). This agreement is obtained thanks to /diss = 1-3. 
In the case of /diss — 1 as in NY04, the models predict 
smaller velocity dispersions than the observations. The rea- 
son /diss — 1-3 is required in this model would be the dif- 
ferences of the SN feedback and the SF time-scale. Dwarf 
galaxies with ctid ^ 20 km s^^ can be obtained only by as- 
sociating with the dynamical response (see Figure|2} as shown 
by NY04, in which the effect of the dynamical response was 
introduced into SA models for the first time. Detailed discus- 
sion on the velocity dispersion of dwarf galaxies can be found 
in the previous subsection and NY04. 

4.5. color-magnitude relations 

Figure ^] shows the V — K color-magnitude (CM) dia- 
gram for cluster elliptical galaxies. Cluster halos that we 
pick out have masses of 6.18, 4.07, 2.99, 2.71, and 2.68 
xlO^*A/0, which are five of the most massive halos in our 
simulation. Different symbols represent elliptical galaxies 
in different clusters. The solid line indicates a CM relation 



18 



NAGASHIMA ET AL. 



M 3 



3.4 
3.2 
M 3 



1 ' 


1 

X □ 




(a) SFB 





_ 
□ 

on 

x'^iffl/; 

'^x O 


^ 1 1 h 


^ ^ 

o<-,^ ■ - ./^ 




(b) WFB 
1 


1 


x** - 
o a" 



I I I , I , I 

-22 -20 -IB -16 



My-5Log(h) 

Fig . 17. — V — K color-magnitude relations of elliptical galaxies in five of 
the most massive clusters in the (a) SFB and (b) WFB models. Different sym- 
bols represent ellipticals in different clusters. The soHd Hne is the observed 
CM relation of the C oma cluster with a correction to the aperture effect by 
IKodama et aTHT99l) 



provided by Kodama et al. (1998), which is corrected for the 
aperture effect on the observational CM relation for the Coma 
cluster obtained by Bower, Lucev & Ellis ( 1992). While the 
models agree with the observation for dwarf ellipticals. My — 
5\og{h) > —18, the predicted CM relations become almost 
flat for giant ellipticals. This trend is similar to previous work 
on SA models whose model parameters are chosen so that 
many othe r observations such as luminosity functions a re re- 
produced c'CoEiliniSQQ^ EiiiiHma&^siH?'^^ The 
mechanism itself to make the slope of the CM relation is 
well understood in terms of the SN feedback. The slope for 
dwarf ellipticals is mainly determined by ahot and the break 
point at which the CM relation becomes flat mainly by Vhnt 
llKauffman n & Charloil998HNagashima & Gouda r200lV 

To improve this situation, the metal enrichment should 
be much more carefully considered because the CM 
relation mainly follows the metallicity-l uminosity rela- 
tion fe odama & Arimoto: ll997t iKauffmaniT & Chariotli ll998t 
iNagashima & Gouda 200 iT In this paper, the chemical yield 
is set at twice the solar value (e.g. Kauffmann & CharloJ 
fT998HSomerville & P rimack 1999). This is required to obtain 
colors of cluster ellipticals similar to observations. However, 
this value is too large compared with yields derived from stan- 
dard IMFs for disk star formation such as Salpeter's, which 
pr ovide almost solar valu es. Moreover, a recent SA analysis 
by'Nagashima & Okamoto (2004) has shown that a Salpeter- 
like IMF successfully reproduces observed stellar distribu- 
tions on an [0/Fe]-[Fe/H] plane. A possible solution is to 



adopt d ifferent IMFs for different modes of star formation. 
iBaugh e t al. (20Q1) and INagashim a et alJ (l2005al) have sug- 
gested Kennicutt's IMF for disk star formation and a top- 
heavy IMF for starbursts. Although this should be worth con- 
sidering, it is beyond the scope of this paper at this stage of 
our modeling. 

5. HIGH REDSHIFT GALAXIES 

5.1. Faint Galaxy Number Counts 

Counting galaxies in a unit are a is a task fundamental to 
observations in astrophysics (e.g.. lYoshii & Takahara 198&) . 
It also provides impo rtant constraints on galaxy for mation 
theories as shown in 'Nagashima et alJ i2001[ |2Qn3)- Fig- 
ure shows faint galaxy number counts in the optical 
BV Ri' z' -h ands. Symbols i ndicate observed counts taken 
from K ashikawa et al.l J2004 I). obtained by the Subaru Deep 
Survey (SDS) using the Suprime-Cam with a substantially 
large field-of-view, 34' x 27', so the data are very reliable in 
a statistical sense, compared with observations by other tele- 
scopes such as the Hubble Space Telescope (HST). In this fig- 
ure, we add offsets to magnitudes depending on pass-bands 
to avoid confusion. Note that these data are observationally 
raw, which means that no correction to selection effects is 
considered and that the isophotal AB magnitude is used. The 
solid and dashed lines represent the SFB and WFB models, 
respectively. The thick and thin lines denote models with 
and without selection effects arising from cosmological dim- 
ming of surface brightness jYoshiafT99l . To estimate the 
selection effects, we need detection limits of surface bright- 
ness dependent on observations and on seeing of observations 
which determines the FWHM of images. The details of the 
selection effects a re sho wn inLTotani & Yoshii (2000), and in 
INagashima et alJ ( 1200 11) within the framework of SA mod- 
els. In these optical pass-bands, as shown in Figure [m both 
models agree incredibly well with the observed galaxy counts 
when the selection effects are included. We have confirmed 
that these models also reproduce galaxy counts well in the 
Hubble Deep Field (HDF), as shown for the Mitaka model in 
NY04. Note that models of dust optical depth for the SFB 
and WFB models are different at high redshift as shown in 
ea.( l38t . We adopt p ~ 1 for the SFB model and p = for 
the WFB model, which means that the dust extinction at high 
redshift in the WFB model is stronger than in the SFB model. 
If we adopted p — Q for the SFB model, it would underesti- 
mate the number of faint galaxies because of the strong dust 
extinction. 

Figure ^] shows galaxy counts in the near-infrared K'- 
band. Open ckcles associated with error bars indicate ob- 
served counts in the Subaru Deep Field (SDF) taken from 
iMaihara etai] (l200l . We should warn readers that the field- 
of-view for the near-infrared SDF is only 2' x 2', which is 
less than 1% of the area for the optical SDF. Other sym- 
bols ar e ground-based observed count s fromlGlazebrook et all 
(Il994t open s t ars), Minezak i et alJ ( 119981 fi lled squares), 
iMcLeod et alJ J 19951 filled triangles), and iGardner et alJ 
j 19961 filled circles). The plotted data are also raw counts. 
The lines denote the same meanings as those in Figure 
ITSl For A''-band counts, moreover, the completeness aris- 
ing from noise and statistical fluctu ations is also taken into 
account. The details are shown in iTotani et"al] ( 120011) and 
Nagashima et al. (2002). Different from the optical counts 
shown in Figure the agreement with the observations is 
worse for both models, which overestimate the number of 
galaxies at K' ^ 23. The SFB model provides somewhat 
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Fig . 18. — Galaxy number counts in the BVRi' z' bands from the left to the right. To avoid any confusion, offsets of -6, -3, +3, and +6 magnitudes are added 
for B, V, i' , and z' bands, respectively. The solid and dashed lines represent the SFB and WFB models, respectively. The thick and thin lines denote models 
with and without the selection effects, respectively, arising from the same detection threshold of suiface brightness of galaxies as employed in the SDS. Symbols 
indicate the observed SDS counts. 
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Fig. 19. — Galaxy number counts in the K' band. Same as Figure fTslbut 
for the near-infrared counts in the SDF. Ope n circles associ ated with error 
bars indicate the observed SDF counts from IMaihara et all f2Q01i.l . Other 
symbols are ground-based observed counts as indicated in the figure. 



better agreement than the WFB model. This is because, in the 
WFB model, more galaxies are formed at high redshift than 
in the SFB model, due to the weakness of the SN feedback. 
Note that the K' band magnitude mainly reflects the stellar 
mass of galaxies, while optical pass-bands magnitudes reflect 
the SFR. Because the K' band magnitude is not significantly 
affected by dust extinction compared with optical bands, one 



might consider that a possible way for the WFB model to re- 
cover the observed near-infrared counts should be stronger se- 
lection effects. However, such a manipulation also affects the 
optical galaxy counts, which are already reproduced. Thus, 
in the range of current modeling, we prefer the SFB model to 
the WFB, from the K' band galaxy counts. 

Finally, we should recall that the above results are due to 
the choice of the redshift dependences of the size and optical 
depth, p and 7. In order to obtain the reasonable agreement 
shown here, it is necessary to introduce these parameters into 
simple models of estimating the disk size and the dust optical 
depth. Although it is quite difficult at this stage to derive the 
values physically, they will provide important keys to under- 
standing of the physics of galaxy formation at high redshift. 

5.2. Redshift distributions 

Figure |20l shows /814 band photometric redshift distribu- 
tion for the HDF galaxies. From the top to the bottom, galax- 
ies for the bright-magnitude bin to the faint-magnitude bin 
are shown. The solid and dashed lines represent the SFB 
and WFB models, respectively. The thick and thin lines de- 
note models with and without the selection effects, which are 
estimated from the same detection criteria as those for the 
HDF. The histograms are the observed pho tometric redshift 
distribution taken from lFurusawa et all ([2000). In the bright- 
est bin (panel a), both models agree well with the observa- 
tion. The SFB model also agrees with the observation over 
all bins. However, as magnitudes of galaxies go fainter, the 
deviation from the observation for the WFB model increases. 
Although the peak of the distribution of galaxies above the 
surface brightness threshold for the HDF is similar to the ob- 
servation, too many galaxies are detected in the WFB model 
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Fig. 20.— Redshift distribution of tlie HDF galaxies for (a) 22 < /814 < 
24, (b) 24 < /gi4 < 26, and (c) 26 < /814 < 28. The soUd and dashed 
lines represent the SFB and WFB models, respectively. The thick and thin 
lines denote models with and without the selection effects, respectively. The 
histogram in each panel is the observed photometric redshift distribution by 
IFurusawa et aL^2000.) . 
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Fig. 21. — Normalized cumulative redshift distribution for the SDF galax- 
ies for (a) 19 < K' < 21, (b) 21 < K' < 22, (c) 22 < K' < 23, and 
(d) 23 < K' < 24. The solid and dashed curves represent the SFB and 
WFB models, respectively. The histograms shown by the solid lines indicate 
the observed redshift distribution. Those indicated by the thin lines show 
the ±3(7 deviated counts estimated by Monte Ca rlo realizations wh en photo- 
metric redshift errors are taken into account Isee lKashikawa et^I] O003t) for 
details]. 



to be consistent with the observation. The weak SN feedback 
is obviously in part responsible for this deviation, because it 
cannot suppress the formation of dwarf galaxies at low red- 
shift. Another important reason is the selection effects. As 
shown by the thin lines, the number of galaxies including un- 
detectable ones is comparable (panel b) or sufficiently large 
(panel c) compared with the observation. However, since 
the selection effects are quite strong, only a small fraction of 
galaxies at high redshift can satisfy the detection criteria for 
the HDF. Therefore, the number of detectable galaxies for the 
WFB model is smaller than that for the SFB model at high 
redshift. Of course, it is possible to increase the number of 
detectable galaxies at high redshift by changing the detection 
criteria and/or by changing the redshift dependence of disk 
sizes. In this case, however, too many faint galaxies would 
emerge to reconcile with the number of galaxies at low red- 
shift and with the observed galaxy number counts. It is evi- 
dent that suppression of galaxy formation is necessary only at 
low redshift in order to match the redshift distribution for the 
WFB model with the observation. Thus, unless such mech- 
anisms are found, the WFB model is inconsistent with the 
observation. 

Figure |^ shows K' band normalized cumulative redshift 
distribution for the SDF galaxies. The solid and dashed lines 
represent the SFB and WFB models, respectively. Only mod- 



els with the selection effects are plotted. The histograms indi- 
cated by the solid lines are the observed redshift distribution 
taken from Kashikawa et al. (2003). Those indicated by the 
thin lines show the ±3(t deviated counts estimated by Monte 
Carlo realizations when photometric redshift errors are taken 
into account. The SFB model agrees significantly well with 
the observation in the top two panels a and b, which shows 
galaxies in the two brightest bins. In panel c, the agreement 
is only marginal. Although in panel d, the deviation from 
the observation is large, we should note quite large errors in 
estimating photometr ic redshift in this magn itude range, as 
shown in Figure 2 of iKashikawa et all J2003 I). On the other 
hand, the difference between the WFB model and the obser- 
vation is large for faint galaxies, K' > 22 (panels c and d). 
The main epoch of galaxy formation in the WFB model is at 
lower redshift than that in the SFB model, as shown in Figure 
|20| 

In IKashikawa et alJ J2003h . the observational data were 
compared with an old version of the Mitaka model given by 
■Nagashima et al, (2£)02). In that model, parameters are simi- 
lar to those of the SFB model, rather than those of the WFB 
model. Nevertheless, the behavior of the old Mitaka model is 
very similar to that of the WFB model. We have confirmed 
that the Mitaka model presented in NY04 also shows a ten- 
dency similar to the WFB model. The reason the redshift dis- 
tribution for the SFB model is improved should be the adop- 
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Fig. 22. — (a) Normalized cumulative redshift distribution for the GOODS 
Ky^g < 22 galaxies, (b) Same as panel a but for the differential redshift 
distribution. The solid and dashed lines represent the SFB and WFB models, 
respectively. The histograms are the observed redshift distribution. Error bars 
in panel b denote 1 a Poisson scatter. To obtain the isophotal magnitude used 
in the GOODS observations, a coiTection of 0.2 mag is applied to the models 
according to iSomerville et alj 120041) . 



tion of merging histories of dark halos directly taken from 
iV-body simulations. This might suggest that using correct 
merging histories is crucial to understanding of high redshift 
galaxies. 

Finally, we compare our models with recent, larger samples 
from the Great Observatories Origins Deep Survey (GOODS) 
data. Figure|22]shows the (a) normalized cumulative redshift 
distribution and (b) differential redshift distribution for the 
GOODS galaxies of KABi= Kycga + 1-85) < 22. The solid 
and dashed curves represent the SFB and WFB models, re- 
spectively. The h istograms are the observ ed GOODS redshift 
distribution from lSomerville et alJJ2004! l Error bars in panel 
b denote 1 cr Poisson scatter. To obtain the isophotal magni- 
tude, according to their paper, we make all galaxies 0.2 mag 
fainter than their total magnitudes. The SFB model seems to 
provide redshift distribution similar to their SA model shown 
in their Figure 1 . Although the SFB model a little overesti- 
mates the number of galaxies at z ^ 1, the predicted number 
of those galaxies at z > 1.5 seems slightly closer to the obser- 
vation than their SA model. In the WFB model, similar to the 
previous figures on redshift distribution, the agreement with 
the observation is worse compared with the SFB model and 
the fraction of low redshift galaxies becomes larger This is 
also because of the weakness of the SN feedback. Thus, we 
can also say that the SFB model is more consistent with the 
GOODS redshift distribution than the WFB model. 
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Fig. 23. — Isophotal area of the SDF galaxies against i^'-magnitude. The 
solid and dashed lines associated with en'or bars of 1 cr scatter represent the 
SFB and WFB models, respectively. To avoid confusion, small offsets ai'e 
added to the magnitude of the WFB model. The isophoto is calculated from 
the observational condition employed in the SDF survey. The crosses indicate 
the observational data from the SDF survey. 

In short, the SFB model is consistent with the photometric 
redshift distribution for /gu < 28 and K' < 23 and the WFB 
only for /814 < 24 and K' < 22. Currently, observational 
data are still limited in number and in accuracy of estimation 
of redshift, particularly for near-infrared pass-bands. Wider 
and deeper surveys in future with accurate determination of 
redshifts, will provide better constraints on galaxy formation. 

5.3. Angular sizes 

Figure |23l shows the isophotal area of galaxies plotted 
against their JiT'-magnitude. The solid and dashed lines as- 
sociated with error bars of 1 ct scatter represent the SFB and 
WFB models, respectively. To avoid confusion, small offsets 
are added to the magnitude of the WFB model. To obtain the 
isophoto, the same observational condition employed in the 
SDF survey is used. The crosses indicate the observational 
data from the SDF survey. This good agreement with the ob- 
servations means that the models correctly estimate not only 
the size of galaxies, but the selection effects arising from the 
cosmological dimming of surface brightness of galaxies. 

5.4. Colors 

The evolution of colors of galaxies depends on growth of 
stellar mass, stellar ages, and dust extinction. Thus we are 
able to extract a great deal of information from galactic colors. 

Figure |23 shows R ~ K colors of galaxies in the AB 
system against their redshift. Because of comparison with 
the GOODS galaxies, we also adopt the selection criterion 
Kab + 0.2 < 22, where we add 0.2 mag to our models 
to obtain the isophotal magnitude. In panel a, the GOODS 
galaxies, retrieved from Figure 2 of ISomerville et all J2004ft . 
are plotted. In panels b and c, galaxies in the SFB and 
WFB models are plotted, respectively. Among galaxies sat- 
isfying the above criterion, we randomly pick out galaxies 
so as to show almost the same number of galaxies as that 
of the GOODS samples. The solid lines indicate evolution- 
ary tracks of single stellar populations of the solar metallic- 
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Fig. 24. — Color-redshift diagram for the (a) GOODS galaxies, (b) SFB model, and (c) WFB model. The color [R — K)yig is defined in the observer frame. 
The observed data are taken from Figure 2 in Somerville et al. 1 2004). In panels b and c, we randomly pick out galaxies of K^b < 22 — 0.2, where 0.2 means 
the correction to obtain the isophotal magnitude, from the mock catalogs based on the SFB and WFB models in order to be a number similar to the GOODS 
samples. The solid lines indicate evolutionary tracks of single stellar populations of solar metallicity with ages of 13.5, 5.8, 3.2, 1, 0.5, and 0.1 Gyr, from the top 
to the bottom, computed using the model given bv .Kodama & Arimoto. i. 1997i) . These ages are the same as those in Figure 2 of .Somerville et al. i2004) . selected 
for ease of comparison. 



ity with ages of 13.5, 5.8, 3.2, 1, 0.5, and 0.1 Gyr, from 
the top to the bottom, estimated using the model given by 
[Kodama & Arimoto ( 199 7i). The se ages are the same as those 
in Figure 2 of S omerville et al.l (|2Q04) so that we can easily 
compare our results with theirs. As can be seen, both the SFB 
and WFB models are in good agreement with the GOODS 
survey. Somerville et al, (.2Q0 4) has concluded that their SA 
model broadly agrees with the GOODS survey, while a small 
discrepancy remains insofar as their SA model does not pro- 
duce enough galaxies with very red colors, (i? — K)ab ^ 4, 
at z ~ 0.8—1.2. Our models seem to be able to pro- 
duce more red galaxies, but slight discrepancies also remain. 
As shown below, one of these is the lack of galaxies with 
{R — K)ab ^ 2 at z ^ 2.5, and another is the existence 
of too many galaxies above the 13.5 Gyr single stellar pop- 
ulation line. A key to understanding the evolution of colors 
should be dust extinction. Our SA model has two types of 
dust extinction associated with the two modes of star forma- 
tion, quiescent disk star formation and starbursts. By switch- 
ing on and off the dust extinction, we investigate the effects 
of dust on colors. 

Here we examine the effe cts of dust extinction for star- 
bursts. As explained in ^ ^2.7\ the SF time-scale for starbursts, 
Tburst, is assumed to be equal to the dynamical time-scale of 
the spheroidal component, <dyn = ''e/Vb. Firstly we change 
this relation to rbmst — O.Sidyn, that is, we take half the dy- 
namical time-scale as the SF time-scale for starbursts. In pan- 
els a and b of Figure |25] we show the same as Figure |23 
but for the above short starburst time-scale for the SFB and 
WFB models, respectively. It is evident that, in both mod- 
els, many galaxies with {R — K)ab ^ 2 at z ^ 2 emerge 
and therefore the agreement with the observation is much im- 
proved, though not perfect. The reason for this improvement 
is explained as follows: The time-scale of starbursts in our SA 
model is distributed in a broad range comparable to the time- 
step, about 0.3-0.5 Gyr, at low redshift. Right after starbursts, 
the age effect on the luminosity, even for such a short time- 
scale, is significant, so we put the epoch of starbursts within 
the time-step randomly, depending on the merger time-scale. 
Therefore, for some galaxies, starbursts happen just before 
the output redshift. Such galaxies just after starbursts are ex- 



pected to be intrinsically very blue. When the SF time-scale 
for starbursts becomes shorter, the amount of dust estimated 
from the amount of cold gas and metals becomes smaller at 
the output redshift. Thus, by shortening the SF time-scale for 
starbursts, several blue galaxies emerge. 

This can be clearly seen if we do not include the dust ex- 
tinction for starbursts at all. In the center panels c and d of 
Figure |25] we show the same models as the SFB and WFB 
models but without dust extinction during starburst phases 
('bare' starbursts). At high redshift z > 1, very blue galax- 
ies with {R ~ K)ab < emerge in both the SFB and WFB 
models. These galaxies are those just after starbursts. The 
GOODS data at such a high redshift suggest that clearly star- 
bursts should not be bare, but our dust model for starbursts 
might be too strong. We have also tried the slab dust model 
instead of the screen dust model, but found that it does not 
work well. It should also be noted that dust extinction can re- 
move galaxies from this figure because galaxies do not satisfy 
the selection criterion, Kab + 0.2 < 22. More sophisticated 
modeling of dusty starbursts would be required in future. 

Next, we discuss the extinction of disk luminosities. In the 
right panels e and /of Figure |25j we show the same models 
as the SFB and WFB models, but without any dust extinc- 
tion. By this manipulation, galaxies above the line of sin- 
gle stellar populations of 13.5 Gyr disappear, similar to the 
GOODS data. This might suggest that our dust model for 
disks is stronger than the observed galaxies. Part of the rea- 
son should be the higher value of the chemical yield we use, 
because the optical depth is proportional to the metallicity of 
cold gas, as mentioned in ii il4.5l Presumably, careful choice 
of IMFs, which determine the chemical yield, depending on 
the SF mode will improve the color distribution. 

It should be noted that the changing of the optical depth for 
disks may affect various observables. For example, weaken- 
ing the dust extinction in a simple way would cause overes- 
timation of faint galaxy counts in optical pass-bands. In or- 
der to obtain full agreement after changing the optical depth, 
therefore, a new set of model parameters will be required. In- 
stead, if the extinction only for strongly reddened galaxies is 
weakened, it might not significantly affect model parameters. 
The change of 7 might also recover the agreement with the ob- 
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Fig. 25. — Color-redshift diagram for the (a) SFB and (b) WFB models with Tburst = 0.5tdyn> the (c) SFB and (d) WFB m odels with no dust extinction 
during starbursts, and the (e) SFB and (f) WFB models with no dust extinction at all. The solid lines are the same as those in Figure l24l 



servational galaxy counts. In all cases, the dust model adopted 
in this paper is very simple and should be refined in future. By 
contrast, the dust model for starbursts is not the case because 
the number fraction of bursting galaxies is smaller than that of 
normal galaxies. As shown in the next subsection, however, it 
might affect the number of high redshift galaxies at z > 3, at 
which a substantial fraction of stars is formed via starbursts, 
while the contribution to faint galaxy counts is small. 

5.5. Cosmic Star Formation Histories 

Recently the cosmic SF history, which is a plot of SFRs in a 
comoving volume against redshift, has been recognized as one 
of the fundamental quantities characterizing galaxy formation 
(|Madau et al. 1996), in spite of large uncertainties in its esti- 
mation caused by, for example, dust obscuration and shape of 
the luminosity function below the detection limit. Figure l26l 
shows the cosmic SFR as a function of redshift. The solid 
and dashed lines represent the SFB and WFB models, respec- 
tively. For comparison, the Mitaka CSF model, in which the 
SF time-scale is independent of redshift similar to the SFB 
and WFB models, is also plotted by the dotted line. Symbols 
are the observational data from surveys indicated in the fig- 
ure. Note that the observational values are corrected, taking 
into account dust extinction, and also that there are still large 
uncertainties originating in undetected faint gala xies at high 
redshift. As shown in Figure 21 of O uchiet al.1 12004), ob- 
servationally obtained cosmic SFRs strongly depend on the 
assumed luminosity function, which is guessed from the data 
for detected bright g alaxies. Th e upper limits for values at 
z > 3 estimated by Ouc hi et aP (|2004) are about a factor of 
two larger than plotted values. These are obtained by inte- 
grating luminosity functions down to the luminosity L ~ Q. 
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Fig. 26. — Cosmic SF histories. The solid and dashed lines represent the 
SFB and WFB models, respectively. For reference, the Mitaka CSF model 
in NY04 is also plotted by the dotted line. Symbols with error bars are the 
observed cosmic SFRs as indicated in the figure. 



Taking into account such uncertainties which are not included 
in plotted error bars, the SFB and WFB models are in good 
agreement with the observations. 

As we have already shown in NY04, the cosmic SFR 
strongly depends on the model of the SF time-scale. If the 
SF time-scale becomes shorter toward higher redshift, for ex- 
ample, proportional to the dynamical time-scale of disks or 
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Fig. 27. — Cosmic SF histories. The solid and short dashed lines represent 
the SFB and WFB models, respectively. The dotted and long dashed lines 
represent the Mitaka CSF and DSF models, respectively, in NY04. The thick 
lines indicate the total SFR, and the thin lines the SFR only for starbursts. 



halos, the peak of the cosmic SFR moves to a higher redshift. 
Because the SFB and WFB models have the same SF model, 
which is independent of redshift, their predicted cosmic SFRs 
are similar The Mitaka CSF model, however, shows different 
cosmic SFRs, in spite of having the same redshift dependence 
as the SF time-scale. This would suggest that it is very impor- 
tant for high redshift galaxies to correctl y est imate merging 
histories of dark halos, as pointed out in 8 ^5.21 

Figurel^shows not only the total cosmic SFR (thick lines) 
but also the SFR only for starbursts (thin lines). The solid and 
short dashed lines represent the SFB and WFB models, re- 
spectively. The dotted and long dashed lines represent the Mi- 
taka CSF and DSF models, respectively. The former has a SF 
time-scale independent of redshift, and the latter a SF time- 
scale proportional to the dynamical time-scale of dark halos, 
which becomes shorter toward higher redshift. Although the 
total cosmic SFRs of the SFB and WFB models are similar 
to the Mitaka DSF model rather than the Mitaka CSF model, 
SFRs for starbursts are quite different. The SFB and WFB 
models predict a very large contribution of starbursts to the 
total SFR up to 15-20% at z - 2 and 25-30% at z - 4, 
while the Mitaka CSF model predicts a contribution of star- 
bursts about 10% at z ~ 2 and 20% at z ~ 4 and the Mi- 
taka DSF model at most 2% at z ~ 1. If the IMF for star- 
bursts is top-heavy, as suggested bv lBaueh et alJ J2005 ) and 
iNagashima et al. (2005a), this significant fraction would af- 
fect many aspects, such as the chemical enrichment and the 
cosmic rates of supernovae and 7-ray bursts. 

6. SUMMARY 

We have developed a semi-analytic (SA) model combined 
with high resolution A^-body simulations, which we call the 
Numerical Galaxy Catalog (i^GC). The minimum mass of 
dark halos is down to 3.O4xlO^Af0, which corresponds to 
Vcirc ~ 20 km s~i at z = and 40 km s^^ at z ~ 4. These 
are enough to resolve the effective Jeans mass, or the filter- 
ing mass, after the cosmic reionization. To examine the ef- 
fects of the resolution on galaxy formation, we also used an- 
other simulation with a box twice the size, in which the min- 



imum mass of dark halos is 2.43 x IO^^A/q corresponding 
to Kirc ^ 40 km s~^ and 80 km s^^. Although this res- 
olution is not significantly worse among the currently avail- 
able cosmological A^-body simulations, we have found that 
it alters luminosity functions, especially the faint-end slope. 
Because SN feedback-related parameters (cthot and Vhot) are 
determined so as to match luminosity functions, worse res- 
olutions should affect values of parameters. We would like 
to stress that the resolution of our SA model is much better 
than that of currently available data provided by SA models 
combined with A-body simulations in the public domain. For 
example, the particle ma ss of dark matter in the GIF project 
iKauffmann et aLlllQQ Qa*) is about 65 times larger than ours, 
and that of the GallCS project (Hatton et al. 2003) is about 27 
times larger. These are even worse than our 140 h^^ Mpc box 
simulation. 

Within observational uncertainties of the shape of luminos- 
ity functions, we have adopted two parameter sets. One is the 
strong SN feedback (SFB) model and another is the weak SN 
feedback (WFB) model. In this paper, we have focused on 
comparison of these models with photometric, structural, and 
kinematical properties of observed galaxies. Both models 
reproduce many observations well, and we slightly prefer the 
SFB model to the WFB model. Below, we summarize the 
results. 

1 . We have explicitly shown that both the SFB and WFB 
models are in reasonable agreement with the following 
observational results: luminosity functions in B and K 
bands, cold gas mass-to-B band luminosity ratios of spiral 
galaxies, H I mass functions, size-magnitude relations for 
spiral and elliptical galaxies, and the Faber-Jackson relation 
at z = 0, and faint galaxy number counts in BVRi'z' bands, 
isophotal area, and cosmic star formation histories at high 
redshift. Note that some of the above observations were 
used to fix the model parameters. For example, the SN 
feedback-related parameters were determined by matching 
the luminosity functions. Table|3lsummarizes observations to 
fix the parameters. 

2. The SFB model is in broad agreement with the observed 
galaxy counts in the near-infrared K' band. Most observed 
redshift distributions also are broadly reproduced by the SFB 
model, except for the faintest galaxies (23 < K' < 24) in the 
Subaru Deep Survey, while the results are much improved 
compared with those of the Mitaka model based on the ana- 
lytic EPS formalism, not on A-body simulations. However, 
since errors caused by uncertainties in estimating photometric 
redshifts are very large for such faint galaxies, this discrep- 
ancy may not be taken so seriously. Future observations with 
accurate determination of redshifts and with a large survey 
area will provide sufficiently better redshift distributions 
for very faint galaxies to make clear whether models are in 
agreement with observations or not. By contrast, the WFB 
model is not consistent with the observed K' band galaxy 
counts in the Subaru Deep Field, while it is in agreement 
with optical (BVRi'z') galaxy counts. This model is also 
inconsistent with redshift distributions, in particular, for faint 
galaxies. 

3. The TuUy-Fisher (TF) relation sUghtly favors the WFB 
model to the SFB model, because the slope of the TF relation 
predicted by the SFB model is steeper than the observed one. 
This contradiction in the TF relation may be improved if we 
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introduce dynamical response on the rotation velocity to gas 
removal induced by the SN feedback, which is only taken 
into account for starbursts making spheroidal systems in the 
current model. We are going to explicitly consider this effect 
on spiral galaxies in future. 

4. The bright end of color-magnitude relations of cluster 
elliptical galaxies in the SFB and WFB models is almost flat 
against magnitudes and also inconsistent with observations. 
Presumably, this will be amended by adopting different 
IMFs for disk star formation and starbursts. Introducing the 
chemical enrichment due to SN e la into the Mitaka model, 
iNagashima & Okamotol i2004h have shown that chemical 
yields consistent with Salpeter's IMF are in agreement 
with observed metallicities and their ratios in the solar 
neighborhood. Moreover, Naaashima et al. (2005a) have 
also shown that metallicities of the intracluster medium 
can be reproduced only when adopting a top-heavy IMF 
for starbursts and Kennicutt's IMF for disk star formation. 
Therefore it is not unreasonable to adopt different IMFs. At 
this stage, however, it is beyond the scope of this paper 

5. Both models predict {R — K)ab colors against redshifts 
broadly consistent with the observational GOODS galaxies. 
Particularly, our models are able to produce red galaxies 
(i? — K)ab 5, which have not been reproduced by an 
SA model of the GOODS group ( Somerville et al. 2004). 
However, slight discrepancies remain. One is lack of galaxies 
with {R - K)ab 1 - 2 at 2; > 2. If we switch off 
the dust extinction for starburst galaxies, several galaxies 
emerge. At the same time, however, very blue galaxies with 
(i? — K)ab ^ 1 also emerge. When a shorter star formation 
time-scale is used for starbursts similar to half the dynamical 
time-scale, the situation can be improved. This means that 
more sophisticated modeling of dusty starbursts should be 
required. Second is that our model produces excessively 
red galaxies at low redshift, z < 1. This is caused by too 



strong a dust extinction for disk stars. This might suggest 
that the chemical yield used here, which is almost twice the 
solar metallicity, is too large, because dust optical depth is 
proportional to the metallicity of cold gas. This high value 
is required to produce red elliptical galaxies to be consistent 
with observed cluster ellipticals. Therefore, similar to color- 
magnitude relations, adopting different IMFs for disk star 
formation and starbursts will help these galaxies become blue. 

This is a first paper in a series of vGC Subsequent pa- 
pers will discuss clustering properties of spatial distribution 
of galaxies ( Yahagi et al. 2005). We also plan to combine the 
quasar formation model given bv'Enoki. Na gashima & Goudal 
( 2003 ) with this model ( Enoki et al. 2005). ^GC wiU be a use- 
ful tool with which to understand galaxy formation in studies 
of various current and future large galaxy surveys. 
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